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PREFACE 


This study under NASA Contract NAS5-21872 was initiated to perform a parametric 
analysis using computer simulation and analysis techniques of the threshold and signal 
distortion effects in FM TV systems. This report is a study of the signal distortion 
effects. The FM threshold study is presented in a separate report. 

Computer simulation is used to model an FM TV transmission system. A T-pulse-and- 
bar test signal is passed through the simulated system and a distorted output signal 
approximating the input test signal is the result of the simulation. The output signal 
can be compared with the input signal for various specified systems. Thus, the dis- 
tortion quality of various systems can be evaluated. 

Asa result of the analysis, the following conclusions can be drawn: 

1. The T-pulse-and-bar test signal is distorted more in a system with preemphasis 
than in one without preemphasis. 

2. In the case of worst distortion a preemphasized system with a half-power RF 
bandwidth equal to twice the peak deviation produces a T-bar overshoot of only 9. 9% - 
less than the 13% limit specified by C. C. I. R. This occurs for a four-pole Chebyshev 
filter having a 0. 1-dB ripple bandwidth of 10. 0 MHz and a half-power bandwidth of 

12. 13 MHz. 

3. An increase in the peak deviation of a system increases signal distortion. 

4. The signal distortion increases as the number of poles of the predetection 
filter increases. 

5. A four-pole Chebyshev predetection filter causes slightly more signal distortion 
than a four-pole Butterworth filter with the same half-power bandwidth. 
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SECTION 1 - INTRODUCTION 


Reception of a frequency-modulated (FM) television signal requires predetection 
filtering to remove noise outside the FM frequency band prior to demodulation. To 
reject as much noise as possible the predetection filter should be as narrow as 
possible without significant distortion of the desired signal. 

The primary objective of this study is to simulate an FM television transmission 
system by computer programming and to evaluate the signal distortion due to 
predetection filtering. The system is modeled mathematically, and a T (T - 100 
nanoseconds) pulse-and-bar test signal is used as the input video signal. The 
output wave resulting from the test signal is calculated and plotted. Then the 
distortion of the signal can be seen and used as a measure of the distortion 
quality of the system. No audio subcarriers are considered in this analysis. The 
simulation is implemented for 64 cases specified by NASA, and the results of those 
cases are presented herein. 

A typical video luminance signal has a format as shown in Figure 1-1, and the 
test signal is shown in Figure 2-3. Portions of the test signal are shown in more 
detail in Figures 4-1 and 4-2. Note that the T-pulse has a width of 100 nanoseconds, 
corresponding to a 5.0-MHz video bandwidth, rather than 125 nanoseconds, correspond- 
ing to a 4. 0-MHz video bandwidth. Therefore, the distortion results will be somewhat 
conservative for application to a 4. 0-MHz system. 

The system being modeled is shown in Figure 1-3. The video test signal enters the 
preemphasis filter which is a 525-line filter specified by C.C.I. R. (Reference 3). 

The model also includes the case for a flat, or unpreemphasized, system in which 
the preemphasis and deemphasis filters are not used. 

The signal goes next to the RF frequency modulator from which the video signal is 
transmitted in an RF form through a noisy RF channel; however, RF additive noise 
is not included in this model. 
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Figure 1-1. Typical Video Luminance Signal 
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Figure 1-2. Block Diagram of the FM Video Transmission System Being Modeled 
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In passing through the predetection filter the signal is distorted somewhat depending 
on the filter characteristics. Narrowing the predetection filter bandwidth cuts out 
noise. On the other hand, narrowing the filter causes signal distortion. Thus, the 
filter is desired to be narrow but not so narrow that significant distortion occurs. 

The predetection filter characteristics are inputs to the computer program for system 
evaluation. 

The RF demodulator recovers an estimate of the modulating signal. It is only an 
estimate because the signal includes distortion. In a flat system, this is the 
estimate of the video test signal. In a preemphasized system the demodulated 
signal must pass through a deemphasis filter to recover the estimate of the video 
test signal. 

The final output signal is compared with the input test signal for various specifications 
of three- and four-pole Butterworth and Chebyshev filters and for various peak 
deviations of the RF carrier. The distortion of the final test signal is used as a 
measure of system quality. 
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SECTION 2 - APPROACH 


2.1 INTRODUCTION 

The mathematical theory required to reclaim a modulating wave in a flat or preemphasized 
FM television system is presented first with block diagrams to illustrate the sequences 
of mathematical operations. Then, after the standard test signal is described, special 
treatment of the test signal to facilitate analysis is explained. Finally, the numerical 
methods used in the analysis are indicated. 

Computer implementation of the mathematics involved in the analysis demands careful 
consideration particularly because convolutions are involved. Convolutions performed 
by numerical methods require lengthy computations in either the frequency domain with 
the use of Fourier transforms or directly in the time domain. Much is written about the 
economy of the fast Fourier transform (FFT) technique, and at first it appears that use 
of the FFT would be more economical than direct implementation of the convolution in the 
time domain. This is true when operations involve the same number of points in both the 
frequency and the time domains. If the number of points required for an operation in the 
time domain is less than that required for the equivalent operation in the frequency 
domain, then the time domain is preferred. To illustrate this point, consider the time 
impulse function, $(t), whose Fourier transform is unity for all values of frequency, 
while <J(t) is nonzero at only one point in the time domain. Some other difficulties with 
the FFT method are integration with respect to time and division of two time functions. 

For these various reasons the time convolution is preferred over the FFT method for 
this application. 

2. 2 A MATHEMATICAL ANALYSIS OF THE VIDEO TEST SIGNAL DISTORTION 

This mathematical analysis involves a method to reclaim the modulating wave of a flat or 
preemphasized FM television system. The case of a preemphasized FM television system 
is first considered, while the case of a flat FM television system can be regarded as a 
special case of the former. 
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Let f(t) denote a test signal which enters the preemphasis filter having impulse response, 
Cjh (t), where h^t) is the response of a passive filter and Cj is the voltage gain required 
to satisfy the preemphasis specification. The modulating wave m(t) after preemphasis is 
given by 

t 

m(t) = c 1 J f<x)h 1 (t - x)dx (2-1) 

— oo 

The modulated carrier which is filtered in the predetection filter is 

s^(t) = cosj<jt + 0(t)j (2-2) 

where u is the angular frequency of the carrier and 
o 

t 

0(t) = k J m(x)dx (2-3) 

— OO 

A constant of proportionality, k, relates the instantaneous frequency deviation to the 
voltage of the modulating wave, m(t). 

The modulated carrier is then passed through the predetection bandpass filter. The 

impulse response, h (t), of the bandpass filter can be found from the specification of the 

z 

filter and can be written for a filter symmetric about its center frequency, , as 

h 2 (t) = 2^(1)008 w Q t (2-4) 

where h f (t) is the impulse response of the equivalent low-pass filter corresponding to 
the bandpass filter. The theory of converting a bandpass to an equivalent lowpass filter 
for analysis purposes is developed in Reference 1. The output, s 9 (t), of the filter is then 
found by convolving the RF signal, s (t), with the impulse response, h 2 (t), as follows: 

t 

s (t) = J cos Lx + 0(x)jh 2 (t - x)dx (2-5) 

-OO 

On substituting for h (t - x) the output becomes 

Z 
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( 2 - 6 ) 


s 2 (t) = J cos[w o x + 0(x)j2h f (t - X) cos [*(t - X)jdx 

— OO 

On expanding the last cosine term above and multiplying one gets 

t 

s ^(t) = cos ujt- J 2 cos w^x cosjo^x + ^( x )Jh^(t - x)dx 


+ sin u> 


t J 2 sin cj q X cos jcu^x + <£(x)jtyt - x)dx 


(2-7) 


A further expansion of the terms of this equation gives 

t 


s 2 (t) = cos <J o t j cos[2<JX + 0(x)j^(t - x)dx 


+ cos ^t j cos ^(xjh^t - x)dx 


+ sin (J 


^t J sinjjiwx + ^(x)jh £ (t - x)dx 


( 2 - 8 ) 


-sin u) t 
o 


t 

j sin <£(x)iyt - x)dx 

— oo 


Note that the first and third integrals are convolutions of an RF signal with the impulse 
response of a lowpass filter. Thus, those terms are zero and s 2 (t) can be written as 

s 9 (t) = A(t) cos cj Q t + B(t) sin u ) Q t < 2-9 ) 


where 



( 2 - 10 ) 


and 



— OO 


( 2 - 11 ) 
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Finally the band-filtered RF carrier can be written as 


s 2 (t) 


VA 2 (t) 


+ B (t) cos 


u t + 
o 


tf'(t)] 


( 2 - 12 ) 


where 


i/>(t) = tan 



(2-13) 


This is a convenient form since the demodulator output is proportional to the derivative 


of ^(t), and 



(2-14) 


The modulating signal, m(t), can be estimated as m(t) by replacing <f>( t) by 0(t) and m(t) 
by m(t) in Equation (2-3). Thus, 

t 

ip( t) = k / m(x)dx (2-15) 


Differentiation of both sides of Equation (2-15) gives ip'(t) - km(t) and the estimate of the 
modulating signal is given by 


m(t) = 


k 


(2-16) 


However, in the preemphasized FM television system m(t) is not yet the final output, 
which is obtainable by passing m(t) through the deemphasis filter having an impulse 
response c h (t), where h (t) is the response of a passive filter and c 3 is the voltage gain 
required to satisfy the deemphasis specification. The final output f(t) is, therefore, 

t 

f(t) = c 3 J m(t)hg(t - x)dx (2-17) 

-OO 


The circumflex accent in f(t) is used to indicate an estimate of f(t). 
m(t) the final output becomes 

A c o * 

f(t) = — f t/,'(x)h 3 (t - x)dx 


On substituting for 

(2-18) 
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A block diagram to summarize the sequence of mathematical operations involved in the 
above analysis is presented in Figure 2-1. 

For the case of a flat FM television system, the modulating wave m(t) is simply the test 
signal, f(t); i. e, , the function, </>{t), in Equation (2-3) is given by 


t 

<£(t) = k J f(x)dx 

-CO 

and the final output, fi(t), is given in accordance with Equation (2-16) as 



(2-19) 


(2-20) 


with k in Equations (2-19) and (2-20) having the same meaning as before. A block diagram 
summarizing the sequence of mathematical operations involved in this flat FM television 
system is presented in Figure 2-2. 

The determinations of the impulse responses of the preemphasis and deemphasis filters are 
described in Appendix A, and those of the equivalent lowpass filters of the two types of 
predetection filters involved are described in Appendix B. 

The determination of the value of k in Equation (2-3) is shown below. To have equal 
maximum frequency deviations on both sides of the center frequency of the predetection 
filter, the zero frequency deviation is assigned at a level 0. 2 volt above the baseline of 
the test signal shown in Figure 2-3. The instantaneous frequency deviation, </>'( t), is given 
by differentiating the expression in Equation (2-3); i. e. , 

(t) - km(t) (2-21) 

Let v be the minimum voltage of the input signal at an instant, t , and let v be the 
1 x z 

maximum voltage of the input signal at another instant, t^. Then, from Equation (2-21) 
the following relations are obtained: 

4 >{ t 1 ) =kv 1 ( 2 - 22 ) 

0 r (t 2 )=kv 2 (2-23) 
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NOTES: 

A: T PULSE 

B: T BAR 

T - 100 ns 


Figure 2-3. T-Pulse-and-Bar Test Signal 
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If D denotes the peak frequency deviation, then 
P 

27TD = 

P 2 


k ( v 2 - V l> 
2 


(2-24) 


for equal maximum frequency deviations from the center frequency of the predetection 
filter. Since the voltage difference, v - v^, is 1 volt, 

k = 47TD (2-25) 

P 


2. 3 THE TEST SIGNAL AND ITS TREATMENT 

The test signal used in the distortion analysis contains a T-puIse and a T-bar, as 

specified in Reference 2. A T-pulse is one cycle of a raised sinusoid having a frequency 

of — . The T-pulse is often referred to as a sine-squared pulse because it can be 
2T 

generated by squaring a sinewave of half the frequency just mentioned. A T-bar is a 
rectangular pulse with half a T-pulse used for the rise and the other half used for the fall. 
The test signal is presented in Figure 2-3. A more detailed view of the pulse and bar 
are shown in Figures 4-1 and 4-2, In any practical situation the T-pulse is sufficiently 
separated from the synchronizing pulse that one can assume settling of any transient 
effect by the beginning of the T-pulse. Similarly one can assume that any transient effect 
preceding the T-bar has settled by the time it starts. The T-pulse is treated as though 
the synchronizing pulse and the T-bar did not exist, and the T-bar is treated as though 
the T-pulse and synchronizing pulse did not exist. Thus, the T-pulse and T-bar are 
considered separately. The methods of decomposing each of these signals to simplify 
the computer simulation are discussed. 

2.3.1 T-Pulse Implementation 

The T-pulse with its dc bias can be represented mathematically as 
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f (t) = 0. 7 Sin 2 w(t - t ) - 0. 2, t sts t 


(2-26) 


= -0. 2 elsewhere 

where a> = (2-27) 

T = 100 nanoseconds 

The duration of the pulse interval (t , t ) is 200 nanoseconds. 

-L Z 

For the case of a flat FM television system, the T-pulse applied to the system is given 
by Equation (2-26) with intervals, (0, t^) and (t^, tg), each of 600-nanosecond duration. 

The signal is considered over the time interval from 0 to t^, This provides an interval 
of observation both before and after the transient pulse. The interval preceding the 
pulse allows time for settling of the system after the signal is applied to the predetection 
filter and before the pulse occurs. The period following the pulse is required to observe the 
delay and distortion on the output signal. The leading interval is required for settling 
of the predetection filter transient which occurs upon application of the signal at time zero. 
The preemphasis and deemphasis filters are tested with a dc voltage level applied for 
a long time prior to the T-pulse. Thus, the resultant dc output signal from these filters 
can be calculated simply as a constant times the input dc voltage without actually using 
the convolution technique. The response to the T-pulse can be superimposed on the 
calculated dc response to arrive at the output waveform. Thus, for the case of a 
preemphasized FM television system, the principle of superposition is applied in order 
to facilitate the analysis by computer. 

The input time function, f(t), is divided up into two components. 


f(t) - f x (t) + f 2 (t) 

where 

f x (t) = 0.7 sln 2 [w(t - tj)], ^ststg 
= 0, elsewhere; 


(2-28) 


(2-29) 
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and 

f (t) = -0. 2, all t (2-30) 

Z 

When f(t) enters the preemphasis filter having an impulse response, c^h^t), the output 
modulating signal, m(t), is given by a convolution operation. 


OO 

m(t) = C 1 J ffxjh^t - x)dx 
-00 


= c i j f 1 (x)h 1 (t - x)dx + c 1 J f 2 (x)h 1 (t - x)dx 


(2-31) 


= c 


1 


/ 


0. 7 sin 


2 



where u is the Heaviside unit function. The upper limit of the integrals is changed to t 
because the filter is causal, and the lower limit of the first integral is changed to 
because the signal, f^(t) , is zero for smaller values of time. The second term on the 
right side of Equation (2-31) is the response of the filter to a dc voltage. If the transient 
effect of applying a dc voltage is assumed to have subsided before the beginning of the 
sine-squared pulse, then the effect of the filter on the dc voltage can be considered as 
merely a constant attenuation. This constant attenuation is specified as a v alue of 10 dB 
in Reference 3. Then, Equation (2-31) can be written as 


m(t) ‘ 3.162278’ t<t l 


(2-32) 


= 1. 479 j 0.7 sin 2 jtu(x - t^j u(-x + t^h^t - x)dx - 3> ^27 8 * t ~ t l 


where the factor, ^227 8’ corresponds 
1. 479, corresponds to a gain of 3. 39 dB. 


to a dc attenuation of 10 dB and the factor. 
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Thus, the modulating signal can be written as the superposition of two signals as follows: 

m(t) = m 1 (t) + m 2 (t) (2-33) 

where 

m^t) = -0. 0632456 for all t 

t (2-34) 

m 2 (t) - l. 479 j 0.7 sin 2 jtj(x - t^uf-x + t^h^t - x)dx, t > 

*1 

- 0 elsewhere 

The first of these signals is the response of the preemphasis filter to the dc portions of 
the input signal, while the second is the response to the sine-squared component. 

Since the frequency modulation and demodulation involved in this analysis are nonlinear 
processes, the principle of superposition cannot be applied as it was with the preemphasis 
filter. The modulating wave m(t) must be applied as it is given in Equation (2-32), but 
the duration of the time intervals can be adjusted appropriately to reduce computation 
time. It was found that a value of 600 nanoseconds for t and a value of 800 nanoseconds 
for the length of the time interval from t to t g were large enough to evaluate the applied 
signal (see Figure 4-131). 

The steps outlined mathematically in Paragraph 2.2 are implemented for finding the 
demodulated output, ^'(t), as an estimate of m(t). The details of this implementation are 
not discussed since no special simplifying technique is used for that portion of the 
calculations. 

Superposition is again applied at the deemphasis filter to simplify the calculations. 

The demodulated output 0 T (t) can be considered to consist of two parts, \p^ andi/^, where 

ip' is the response to the dc portion of the modulating wave, m^(t), and ip 2 is the response 

to m (t). Both m (t) and m (t) are defined by Equation (2-34) . The demodulated output, 

2 12 

t//(t), can be written as 

rp'{ t) = ^(t) + ip' 2 (t) (2-35) 
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consisting of the dc component, ^(t), and the transient component, i/r^t). The input wave 
can be written as 

0(t) = ip^H + (2-36) 

where ^(<*>) is the steady-state dc response having a value of -0, 0632456. 

If c h (t) denotes the impulse response of the deemphasis filter, the response of the 
3 3 

deemphasis filter to the input wave, <// r (t), is given by 

QO 

f(t> = c 3 / fen + V ” x)dx 

-CC L 

t t 

= c 3 J- (0. 0632456)h 3 (t - x)dx + c^J (xjlyt - x)dx (2-37) 

t x 

A A 

- f x (t) + f 2 (t) 

where the upper limit of the integrals has been changed to t because the deemphasis filter 
is causal, and the lower limit of the second integral has been changed to t^ because the 
function, i|f * (t) , is zero for smaller values of time. 

A 

If only the steady-state value, f^t), is considered, it is given by 

£(t) = f (oo) = (-0. 0632456) (3. 162278) (2-38) 

where 3. 162278 is the value of the deemphasis amplification factor, c . 

3 

A 

The final output, f(t) , of the system is given by 

t 

f(t) = -0.0632456(3. 162278) + 3. 162278 I ip' (x)h (t - x)dx 

J A o 

1 (2-39) 

t 

= -0. 2 + 3. 162278 j 

l l 
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From the previous derivations, it is seen that the response of the preemphasis and 
deemphasis filters given by convolution operations are broken up into two separate con- 
volutions, one of which can be evaluated without performing the operation. In this way, 
the required computer time is greatly reduced. 

2.3.2 T-Bar Implementation 

The procedures described for the T-pulse implementation are also applicable to the T-bar 
test signal when it is resolved into the right component parts. The following shows how 
this can be done. 

2. 3. 2. 1 The Rise Portion of the T-Bar 

Let T (t) represent the rise-portion of the T-bar having a mathematical representation 
R 

given by 


T R (t) = -0.2 


-M 

VI 


2 

= 0. 7 sin 

‘|w(t - t^J-0. 2 * 

t.s t<t 
4 5 

(2-40) 

= 0. 5 


t st<t 
5 6 



where t. is the instant at which the half sine-squared pulse starts, 

4 

t is the instant at which the half sine-squared pulse ends, 

5 

the duration of the interval (t., t ) is 100 nanoseconds, 

4 5 

The duration of the interval, (t , t ), is arbitrary as long as it is large enough for the 

5 6 

transient response of the filters to subside. 

To make T (t) suitable for the previously developed procedures of analysis in the case of 
R 

the T-pulse, T (t) is reconstructed in the following way. 

R 

T R (t) = T Rl (t) + T R2 (t) (2-41) 

where 
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t ri« - 0 


t < t A 
4 


2 

= 0.7 sin 

^(t - t 4 )], 

' »4 £t£ *5 

(2-42) 

= 0.7 





and 

T R2 ( t ) = -0. 2, all t (2-43) 

When comparison is made with Equation (2-28), it is clear that T R1 (t) corresponds to 

f (t) while T (t) corresponds to f (t), and the same procedures of analysis can be applied. 
1 ' R2 4 

2. 3. 2. 2 The Fall Portion of the T-Bar 

It is assumed that any transient filtering effects have subsided before reaching the fall 
portion of the T-bar signal. This is reasonable because the pulse is quite long. There- 


fore, the fall portion of the T-bar signal, T (t), can 

be represented by 


T F (t) = 0.5 

t<t ? 


= 0. 5 - 0. 7 sin 2 jw(t - t 7 )J 

• ‘7 St:St S 

(2-44) 

= -0.2 

t>l 8 



where the duration of the interval (t 7 ,tg) * s a lso 100 nanoseconds, and the value of t^ 
is arbitrary. 

T (t) is also reconstructed in the following manner: 

F 

T p (t) = T F1 (t) + T F2 (t) (2-45) 

where 

T F1 (t) = 0 t<t ? 

= -0. 7 sin 2 [(J(t - t 7 )| , t ? <t <t g (2-46) 

= -0. 7 t > t 8 
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and 


T F2 (t) = 0.5, all t 


(2-47) 


Again, with reference to Equation (2-28), it is seen that T (t) corresponds to f^t) and 

T (t) to f (t) so that the same procedures of analysis can be applied. 

F 2 2 

For the case of a flat FM television system, the constant voltage portion of the T-bar 
test signal was shortened to 600 nanoseconds and then the entire modified T-bar signal 
was applied directly to the system. 

2.4 NUMERICAL METHODS USED IN THE COMPUTER IMPLEMENTATION 


The numerical formula used in the analysis for evaluating a definite integral is simply 
the rectangular formula. 


2 N 

f f(x)dx f(t + n4t)At 

V n=l 


(2-48) 


where f(t) is a given continuous time function, 

t and t are the limits of integration, 

1 2 

N is the number of sections into which the interval of integration ( t - t^) is 
partitioned, 

l 2 ' t l 

4t = — — — , the size of an increment of time. 


The numerical formula used for differentiation is a three-point one (Reference 4) as 
follows: 


f! = -rr; (-3f. + 4f. 
l 2At ' l l+l 


f i+ 2> 


(2-49) 


where f , f , and f _ are the values of the given function f(t) at the i th , (i + l) th , and 
i i+1 i+2 

(i + 2) th points of the t-axis, 

th 

f' is the approximate value of the derivative of f(t) at the i point of the t-axis, 
At is the step-size into which the t-axis is divided. 
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The numerical formula for evaluating a given convolution of two causal time functions is 
given below, and its rigorous proof can be found in Reference 5. 

Let t 

R(t) = J f(x)g(t - x)dx (2-50) 

o 

where f(t) and g(t) are two causal functions. 

If the range of integration is divided into n intervals of dt second, then R(n4t) is given by 

n 

R(n4t) = At } f g (2-51) 

' ' mn-m 

m=0 


2-17 



SEC TTON 3 - DESCRIPTION OF COMPUTER PROGRAMS 


3.1 INTRODUCTION 

The computer programs* developed in FORTRAN for the FM distortion analysis, are 
discussed in this section and a listing of the programs is found in Appendix C. 

The programs are not combined into a single program, because some results are 
calculated only once and are stored in data files for repeated use as input data for 
some of the other programs. There are sixteen numbered programs each of which 
performs a particular function as shown in Table 3-1. The final results of each 
program are stored in data files which can be printed out if desired or used 
directly as input data for other programs. Some of the computer programs are 
very similar and will be described as a group. 

The computer programs all use a 1-nanosecond time increment for sampling 
and mathematical operations such as integration and differentiation. 

3.2 COMPUTER PROGRAM 1 

This computer program is used to generate sample values of components of a 
T-bar test signal. Its flowchart is shown in Figure 3-1. The data generated 
by this program are used as input data for Program 7. The program calculates 
samples of the rise and fall components of the T-bar according to Equations 
(2-42) and (2-46), respectively. 

It generates samples at one-nanosecond intervals from the leading edge of the 
sine-squared curve at time, Tl, to some later time, TEND, following the end of 
the sine-squared curve. In this program, Tl is always set to zero, and the end 
of the sine-squared curve is designated T3 with a value of 100 nanoseconds. 

The sample values generated by this program are printed out and stored in two 
data files for later use - one for the rise and one for the fall. 
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3.3 COMPUTER PROGRAM 2 


This computer program is used to generate sample values of an unbiased T-pulse. 
The flow chart, similar to that for the T-bar, is shown in Figure 3-2. The sample 
values are calculated, according to Equation (2-29), at 1-nanosecond intervals from 
the leading edge of the sine-squared curve at time, Tl, to some later time, TEND, 
following the end of the sine-squared pulse. The sample values for times later 
than T3 are zero. In this program, Tl is always set to zero, and the end of the 
sine-squared curve is designated T3 with a value of 200 nanoseconds. The sample 
values are printed out and stored in a data file for later use in Program 8. 

3.4 COMPUTER PROGRAM 3 

This program generates samples of the biased and shortened T-bar test signal. 

The flowchart for this program is shown in Figure 3-3. The data generated by 
this program are used as input data for Program 9. It generates samples at 
1-nanosecond intervals from zero time to some later time, TEND (1800 nano- 
seconds). The samples are calculated for the T-bar signal starting at time 
Tl (600 nanoseconds) and ending at T4 (1400 nanoseconds), and they include 
values of the -0. 2-volt dc bias for 600 nanoseconds preceding the shortened 
T-bar and 400 nanoseconds following it. 

The sample values generated by this program are printed out and stored in a 
data file for later use. 

3. 5 COMPUTER PROGRAM 4 

This program generates samples of the biased T-pulse test signal. The flow- 
chart for this program is shown in Figure 3-4. The data generated by this 
program are used as input data for Program 9. It generates samples at one- 
nanosecond intervals from zero time to some later time, TEND (1400 nano- 
seconds). The samples are calculated for the T-pulse signal starting at 
time T2 (600 nanoseconds) and ending at T3 (800 nanoseconds), and they include 
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values of the -0. 2-volt dc bias for 600 nanoseconds preceding the T-pulse and 600 nano- 
seconds following it. 

3. 6 COMPUTER PROGRAMS 5 and 6 

These two computer programs are used to generate sample values of the preemphasis 
and deemphasis filters, respectively. Their flowcharts are identical and therefore 
only one flowchart is shown for them in Figure 3-5. 

The samples, calculated by Program 5 and used as input data for Program 7, are 
printed out and stored in a data file for later use. Those calculated by Program 6 and 
used as input data for Programs 15 and 16 are also printed out and stored in a separate 
data file. 

3. 7 COMPUTER PROGRAM 7 

This computer program is used to generate sample values of the modulating wave for the 
case of a preemphasized FM television system with a T-bar test signal. Rise and fall 
sections of the T-bar are calculated separately and are stored in two separate files. 

The flowchart for this program is shown in Figure 3-6. The input data are generated 
by Programs 1 and 5. Program 1 supplies the samples of the rise and fall components 
of the T-bar test signal, while Program 5 supplies sample values of the preemphasis 
impulse response. This program performs the convolutions required and adds in the 
dc voltage levels to generate samples of the output waveforms of the preemphasis filter. 
These data are printed out and stored in files for later use in Program 9. 

3. 8 COMPUTER PROGRAM 8 

This computer program is used to generate sample values of the modulating wave for the 
case of a preemphasized FM television system with a T-pulse test signal. The flow- 
chart for this program is shown in Figure 3-7. The input data are generated by 
Programs 2 and 5. Program 2 supplies samples of the T-pulse, while Program 5 
supplies sample values of the preemphasis impulse response. This program performs 
the convolutions required and adds in the dc voltage levels to generate samples of the 
output waveform of the preemphasis filter. These data are printed out and stored in 
files for later use in Program 9. 
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3.9 COMPUTER PROGRAM 9 


This program calculates sample values of cos <p( t) and sin <£(t) for both flat 
and preemphasized FM television systems. The flowchart is shown in Figure 3-8.. 
The input samples of the modulating wave generated by Program 3, 4, 7, or 8, are 
used to calculate samples of <f>( t) in accordance with Equation (2-3) for the pre- 
emphasized case and in accordance with Equation (2-19) for the flat case. After 
integrating to obtain <f>( t) the sample values of cos <f>(t) and sin c/>(t) are calculated and 
stored in two separate files for later use in Program 14, 15, or 16. 

3.10 COMPUTER PROGRAMS 10 THROUGH 13 

All these computer programs are used to generate the sample values of the impulse 
responses of three - and four-pole Butterworth lowpass filters, the three- and four- 
pole Chebyshev lowpass filters. The single flowchart for Programs 10 and 12 is 
shown in Figure 3-9, and the flowchart for Programs 11 and 13 is shown in 
Figure 3-10. The variable, FB, in these programs represents the bandwidth in Hertz 
and should be assigned a value accordingly. In Programs 10 and 11, FB is the half- 
power bandwidth of a Butterworth filter, while in Programs 12 and 13 it represents 
the ripple bandwidth of a Chebyshev filter. The resulting sample values of the 
impulse response are stored in a file for later use in Program 14, 15, or 16. 

3.11 COMPUTER PROGRAM 14 

This computer program performs the final mathematical operations required to 
obtain sample values of the final output of a flat FM television system. The flow- 
chart is shown in Figure 3-11. There are three input data files for this program, 
and they contain sample values of cos <£{t), sin <f>(t) and h (t). Values of cos (j>(t) 
and sin <p(t) are generated by Computer Program 9 and values of h|(t) are generated 
by Programs 10, 11, 12, or 13, depending on the predetection filter specifications. 
The peak deviation, DP, must be specified within the program. The program 
functions as follows: 


3-4 



Two convolution operations produce the time functions, A(t) and B(t), defined by 
Equations (2-10) and (2-11). A(t) and B(t) are obtained by convolving cos </>(t) 
and sin <£(t) with the impulse response h|(t) of an equivalent lowpass predetection 
filter. The sample values of cos <jb(t), sin <j>(t), and h^(t) are read from files 
into the vectors FLCOS, FLSIN, and HBF, respectively. The resulting sample 
values of A(t) and B(t) are retained in the vectors A and B, respectively. 

The time function i//(t), which is the arctangent of B(t)/A(t), is then produced and the 
sample values of ^/(t) are retained in the vector PHI. When A(t) takes a value of 
zero, the corresponding calculation of the arctangent of B(t)/A(t) is skipped. This 
fact is reported by printing out the message "A=0.0. " 

Finally, the differentiation of t|»(t) is performed and the sample values of the 
resulting function, ip'(t), are then divided by K. Then the resulting sample values 
of the final output are retained in a file and are also printed out. 

3.12 COMPUTER PROGRAMS 15 AND 16 

Programs 15 and 16 are used to calculate sample values of the final output waveform 
for a preemphasized FM television system. Program 15 calculates the system 
response to the T-pulse or to the rise portion of the T-bar test signal. Program 16 
calculates the system response to the fall portion of the T-bar. These two programs 
are identical except for the values of two constants; therefore, only one flowchart,, 
shown in Figure 3-12, is used to describe the programs. There are four input data 
files for this program, and they contain sample values of cos </>(t), sin 4>(t), h^(t), 
and h 3 ( t) . The values of cos <£(t) and sin <f>(t) are generated by Computer Program 9 
and values of hg(t) are generated by Program 10, 11, 12, or 13, depending on the 
predetection filter specification. The peak deviation must be specified within 
each program. The value M in Program 15 is specified as 1200 for T-bar calcu- 
lations and as 1400 for T-pulse calculations. The value of M does not change from 
1200 in Program 16. Each program functions as follows: 
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Two convolution operations produce the time function, A(t) and B(t), defined by 
Equations (2-10) and (2-11). A(t) and B(t) are obtained by convolving cos <f>( t) 
and sin ^»(t) with the impulse response, (t), of an equivalent lowpass predetection 
filter. The sample values of cos <£(t), sin <£(t), and h^(t) are read from files into 
the vectors FLOCOS, FLSEN, and HBF, respectively. The resulting sample values 
of A(t) and B(t) are retained in the vectors A and B, respectively. 

The time function, i/»{t), which is the arctangent of B(t)/A(t) is then produced and the 
sample values of t//(t) are retained in the vector PHI. When A(t) takes a value of 
zero, the corresponding calculation of arctangent of B(t)/A(t) is skipped. This 
fact is reported by printing out the message "A=0. 0". 

The differentiation of i>(t) is then performed and the sample values of the resulting 
function, «//(t), are divided by K and retained in the vector PHID. Then all these sample 
values are examined for computational errors. The errors, occurring in pairs, are 
corrected by interpolation. 

After corrections are made to the components of the vector PHID, this vector is 
convolved with the vector HDE, which contains the sample values of the impulse 
response of the deemphasis filter, to give the final output of the preemphasized 
FM television system. 


3-6 



Table 3-1. Functions of Computer Programs 


PROGRAM 

NUMBER 


1 


2 


3 

4 

5 

6 
7 


8 


9 

10 


PROGRAM 

FUNCTION 


Generates the sample values of a T-bar sig 
nal for a preemphasized FM television sys- 
tem. 

Generates the sample values of a T-pulse 
signal for a preemphasized FM television 
system. 

Generates the sample values of a T-bar sig- 
nal for a flat FM television system. 

Generates the sample values of a T-pulse 
signal for a flat FM television system. 

Generates the sample values of the impulse 
response of the preemphasis filter. 

Generates the sample values of the impulse 
response of the deemphasis filter. 

Calculates the sample values of the modu- 
lating wave corresponding to a T-bar sig- 
nal in a preemphasized FM television sys- 
tem. 

Calculates the sample values of the modu- 
lating wave corresponding to a T-pulse test 
signal in a preemphasized FM television sys- 
tem. 

Calculates the sample values of the time 
function cos cf>( t) and sin <f»(t) . 

Generates the sample values of the impulse 
response of a symmetrical 3 -pole Butter- 
worth filter. 
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Table 3-1. Functions of Computer Programs 


PROGRAM 

NUMBER 

11 

12 

13 

14 

15 

16 


PROGRAM 

FUNCTION 

Generates the sample values of the impulse 
response of a symmetrical 4-pole Butter- 
worth filter. 

Generates the sample values of the impulse 
response of a symmetrical 3-pole Chebyshev 
filter. 

Generates the sample values of the impulse 
response of a symmetrical 4-pole Chebyshev 
filter. 

Calculates the sample values of the final out- 
put of a flat FM television system. 

Calculates the sample values of the final out- 
put of a preemphasized FM television system 
when the input is either a T-pulse or the rise 
portion of the T-bar signal. 

Calculates the sample values of the final 
output of a preemphasized FM television 
system when the input is the fall portion of 
the T-bar signal. 
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Figure 3-1. Flowchart for Computer Program 1 
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Figure 3-3. Flowchart for Computer Program 3 


3-11 
















Figure 3-4. Flowchart for Computer Program 4 
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Figure 3-5. Flowchart for Computer Programs 5 and 6 
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Figure 3-6. 


Flowchart for Computer Program 7 
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Figure 3-7. Flowchart for Computer Program 8 
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Figure 3-8. Flowchart for Computer Program 9 
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Figure 3-9. Flowchart for Computer Programs 10 and 12 
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Figure 3-10. Flowchart for Computer Programs 11 and 13 
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Figure 3-11. Flowchart for Computer Program 14 
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SECTION 4 - ANALYSIS OF COMPUTER RESULTS 


4.1 INTRODUCTION 

The computer programs described in Section 3 and listed in Appendix B are developed to 
calculate responses of FM television systems to a T-pulse-and-bar test signal. The 
responses are calculated for 64 different combinations of system parameters. These 
responses are used as a measure of the distortion quality of the systems being evaluated. 

The response of a system can be compared with the input test signal and with the responses 
of other systems. 

Table 4-1 is convenient for finding the figure showing the calculated system response for 
particular system parameters. 

4. 2 INPUT TEST SIGNAL 

The T-pulse-and-bar test signal is shown in Figure 2-3 in terms of the time relation of 
the pulse and bar to a single television line scan. Figure 4-1 shows the T-pulse in 
detail, while Figure 4-2 shows the T-bar in detail. The width of the pulse at half ampli- 
tude is T, having a value of 100 nanoseconds. The entire length of the pulse is 200 nano- 
seconds. The height of the pulse is 0. 7 volt extending from -0. 2 volt to 0. 5 volt. 

The T-bar has a width of 25. 0 microseconds at half amplitude and a total length of 25.2 
microseconds. Its height is also 0. 7 volt extending from -0. 2 volt to 0. 5 volt. Note that 
a break in the graph of Figure 4—2 is used to indicate that the T— bar continues at the same 
level over a long period of time from one side of the break to the other. This technique 
is used to show detail of the T-bar on a single page. 

4.3 PRESENTATION OF THE RESULTS 

The calculated system responses to the T-pulse-and-bar are shown in Figures 4-3 through 
4-130. Table 4-1 is an index for these figures. In all cases the T-bar response for a partic- 
ular specified system is found in the figure immediately following the T-pulse response 
for that same system. This is clearly presented in Table 4-1. For example, the T-pulse 
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response of a preemphasized FM system with a three-pole Butterworth predetection 
filter having a half-power bandwidth of 10 MHz and with a 5-MHz peak deviation is 
found in Figure 4-67. The T-bar response for the same system is found in Figure 4-68. 

The RF bandwidths of the predetection filters are specified in two ways. Half-power 
bandwidths are specified for the Butterworth filters and ripple bandwidths are specified 
for the Chebyshev filters. Table 4-2 shows the half-power bandwidths corresponding to 
the ripple bandwidths for all the cases analyzed. 

4.4 DISCUSSION OF THE RESULTS 

The distortion of the test signal in passing through the system can be evaluated as a 
function of the system parameters. The effect of the RF predetection filter type, its 
number of poles, its bandwidth, the peak deviation, and preemphasis are discussed. 

4.4.1 Effects of Predetection Filter Types 

A comparison is made of the distortion effects resulting from the use of Butterworth 
and Chebyshev filters with the same number of poles and the same specified bandwidths. 
Comparison of the cases exhibiting a rather large distortion is first made. Consider 
the cases for the three-pole, 10.0-MHz filters and 5.0-MHz peak deviation. Table 4-1 
shows that Figures 4-3 and 4-4 give the test signal response for a flat system utilizing 
a Butterworth filter and that Figures 4-19 and 4-20 give the corresponding test signal 
response for a flat system utilizing a Chebyshev filter. The undershoot of the T-pulse 
in these cases is 3. 1%(0. 022 volt) for the Butterworth and 2.3% {0.016 volt) for the 
Chebyshev cases. These undershoot values compare, respectively, with 6.9% (0. 048 volt) 
and 2. 1% (0.015 volt) for a preemphasized system. These latter numbers are derived 

from Figures 4-67 and 4-83. The T-bar responses also show more distortion for the 
Butterworth than for the Chebyshev cases. If the same comparison is made for four- 
pole filters, the Chebyshev filter again looks better. On the basis of this comparison of 
the use of a Butterworth filter having a given half-power bandwidth with a Chebyshev filter 
having a ripple bandwidth of the same value, the Butterworth filter produces more distortion. 
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Since this comparison may not be fair to the Butterworth filter a similar comparison is 
made comparing Butterworth and Chebyshev filters having the same half-power band- 
widths as well as the same number of poles. Table 4-2 shows that a four-pole Chebyshev 
filter having a 10. 0-MHz ripple bandwidth has a half-power bandwidth of 12. 13 MHz; 
therefore, the four-pole Chebyshev filter with a 10. 0-MHz ripple Bandwidth can be 
compared with the four-pole Butterworth filter with a 12. 0-MHz half-power bandwidth. 
The comparison for the T-pulse in a flat system can be made by comparing Figures 4-43 
and 4-57. These cases utilize a peak deviation of 6.0 MHz. The undershoot is 3.9% 

(0. 028 volt) for the Butterworth case compared with 4. 3% (0. 030 volt) for the Chebyshev. 
This difference is rather small. On the fall side of the T-pulse the overshoot is 2.9% 

(0. 020 volt) for the Butterworth filter compared with 4.3% (0. 030 volt) for the Chebyshev. 
Thus, it appears that the Chebyshev filter still produces slightly more distortion than the 
corresponding Butterworth filter. When the T-bar signals are compared in Figures 4-44 
and 4-58, one sees that the Butterworth filter produces an overshoot of 2. 7% (0. 019 volt) 
on the fall compared with 2. 3% (0. 016 volt) for the Chebyshev. This would indicate that 
the Butterworth filter produces slightly more distortion than the Chebyshev. Thus, there 
appears to be no significant difference in the performance of the comparable Butterworth 
and Chebyshev filters for a flat system. 

A similar comparison for a preemphasized system shows a more significant difference 
in performance of the two types of filters. A comparison of T-pulse distortions is made 
by observing Figures 4-107 and 4-121. The undershoot is 5. 6% (0. 039 volt) for the 
Butterworth case compared with 6.0% (0. 042 volt) for the Chebyshev, and on the fall 
side of the T-pulse the overshoot is 4.4% (0.031 volt) for the Butterworth compared with 
9. 9% (0. 069 volt), the largest pf all overshoots, for the Chebyshev. That the Chebyshev 
filter produces more distortion than the Butterworth in this case is strengthened by 
comparing the overshoots at the fall ends of the T-bar output signals. Figures 4-108 
and 4-122 show overshoots of 1. 4% (0. 010 volt) and 3. 6% (0. 025 volt), respectively, 
for the Butterworth and Chebyshev filters. Thus, the Chebyshev overshoot at the end 
of each signal is approximately twice as large as for the corresponding Butterworth 
overshoot. 
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For broadband cases the distortion from one filter to another is small enough that any 
differences in the distortions are also small. In such cases the differences may be 
within the accuracy of the computer estimation of the output signal. 

In conclusion, if preemphasis is used, it appears that the Chebyshev filter produces more 
distortion than a Butterworth filter with the same number of poles and the same half-power 
bandwidth and that without preemphasis there is no significant difference. It also appears 
that the Butterworth filter with a certain half-power bandwidth produces significantly more 
distortion than a Chebyshev filter with a ripple bandwidth numerically equal to the half-power 
bandwidth of the Butterworth filter. 

4.4.2 Effects of the Number of Poles of the Predetection Filter 

If all other parameters of the FM system remain fixed, the number of poles of the RF 
predetection filter affect the distortion of the system. In fact, distortion increases with 
the number of poles, and some examples are chosen to illustrate this point. 

First, an example of Butterworth filters is chosen. Figures 4-3 and 4-9 are compared 

for T-pulse distortion. The peak deviation is 5. 0 MHz and the half-power bandwidth is 

10. 0 MHz for this case. The undershoot at the top of the pulse is only 3. 1% (0. 022 volt) for the 

three-pole filter compared with3.7% (0.026 volt) for the four-pole filter. The overshoot at the 

end of the pulse is 4.3% (0. 030 volt) compared with 5. 7% (0.040 volt) with the larger number 

of poles causing the larger distortion. A similar effect is found in a preemphasized system. 

It also shows up in the T-bar responses, which can be compared in Figures 4-4 and 4-10. 

An example chosen for the Chebyshev case is the comparison of Figures 4-99 and 4-105 
for a preemphasized system with a peak deviation of 6. 0 MHz and a ripple bandwidth of 
10.0 MHz. The T-pulse exhibits an undershoot of 7.4% (0.052 volt) at the peak for the three- 
pole filter compared with 8. 6% (0. 060 volt) for the four-pole case. Again a four-pole filter is 
seen to cause a larger distortion than a three-pole filter. 

4. 4. 3 Effects of the Bandwidth of the Predetection Filter 

The bandwidth of the predetection filter definitely affects the distortion of the test 
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signal and this effect is clearly seen in the output waveforms. If all other system para- 
meters are unchanged, an increase in bandwidth results in a decrease in distortion. 

As the bandwidth of the three-pole Butterworth predetection filter increases from 10.0 
MHz to 12. 0 MHz and then to 15. 0 MHz, the T-pulse distortion shown in Figures 4-3, 

4-5, and 4-7 decreases. The T-pulse undershoot at its peak is 3.1% (0. 022 volt) for a 10. 0- 
MHz bandwidth, 2. 6% (0.018 volt) for a 12. 0-MHz bandwidth, and 1. 4% (0, 010 volt) for a 15. 0- 
MHz bandwidth. The overshoot at the end of the pulse decreases from 4. 3% (0. 030 volt) to 
2. 1% (0. 015 volt) and finally to'l. 1% (0. 008 volt) as the bandwidth increases. Similar 
decreases in distortion occur in other examples as the bandwidth of the predetection filter 
increases.. As it increases without other changes, more of the RF spectrum is passed by 
the filter and this results in less distortion of the modulated waveform, which, in the end, 
results in a less distorted output signal. 

4.4.4 Effects of Peak Deviation 

A change in the peak deviation affects the distortion if all other system parameters 
remain unchanged. The figures listed at the top of Table 4-1 can be compared on a 
one-to-one basis with those at the bottom of the table. Those at the top are for a 5. 0- 
MHz peak deviation, and those at the bottom are for a 6. 0-MHz peak deviation. In com- 
paring Figure 4-3 with Figure 4-35 one sees that the T-pulse undershoots by 3.1% (0.022 
volt) at the peak in the first figure and by 4.3% (0. 030 volt) in the latter. However, the 
overshoot at the end of the pulse is 4.3% in both figures. In general, such comparisons 
as this indicate more distortion with more peak deviation. This is to be expected, since 
more of the FM spectrum is rejected by the predetection filter as the peak deviation 
increases. 

4.4.5 Effects of Preemphasis 

If all other system parameters are unchanged, the distortion of the test signal depends 
on whether or not preemphasis is used. It is quite clear from the data that preemphasis 
increases distortion of the test signal. As an example, compare the distorted T-pulses 
of Figures 4-3 (flat) and 4-67 (preemphasized). The undershoot of the pulse peak in the 
flat system is 3. 1% (0. 022 volt) compared with 6. 9% (0. 048 volt) in the preemphasized 
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system, and the overshoot at the end of the pulse is 4.3% (0. 030 volt) in the flat system 
compared with 4. 6% (0. 032 volt) in the preemphasized system. 

This effect can be explained if the modulating signals are compared for the two cases. 

In the flat case the modulating signal is just the T-pulse-and-bar as shown in Figures 4-1 
and 4-2. The modulating voltage ranges from -0. 2 volt to 0. 5 volt and the peak deviations 
correspond to ±0. 5 volt. If a preemphasis filter is added, the modulating waveform 
takes the form shown in Figures 4-131 through 4-133. 

Figure 4-131 shows the response of the preemphasis filter to the T-pulse. Note that 
the modulating voltage reaches a peak of 0. 79 volt corresponding to a peak deviation 
that is 58% larger than that which occurs without preemphasis. Naturally the pre- 
detection filter distorts this peak, since the RF frequency corresponding to the peak is 
farther from the center of the predetection filter. At the end of the pulse the preempha- 
sis filter causes an overshoot to -0. 33 volt compared with -0. 20 volt for the flat case. 
Thus, one could possibly expect more distortion at the end of the pulse in a narrow-band 
system, although the peak deviation is not excessive on that side of the carrier frequency. 

Figure 4-132 shows the response of the preemphasis filter to the rise portion of the 
T-bar. Again, as in the case of the T-pulse, the voltage rises to 0. 79 volt which provides 
a 58% increase in peak deviation over that without preemphasis. Then the voltage 
slowly falls to a lower constant dc level. This results in an increased deviation on one 
side of the carrier frequency and distortion by the predetection filter. 

Figure 4-133 shows the response of the preemphasis filter to the fall portion of the T-bar. 
Here the voltage falls to approximately -0.7 volt, causing an increase of 40% in peak 
deviation over that achieved at the peak of the synchronizing pulse in the absence of pre- 
emphasis. Preemphasis increases the deviation of the T-bar test signal on the low side 
of center frequency by 350%, since the T-bar test signal goes only as low as -0. 2 volt 
in the absence of preemphasis. 

Clearly preemphasis causes increased distortion of a T-pulse-and-bar test signal over 
that of a flat system, and that increased distortion results from the increased deviation 
of the RF carrier and resultant distortion of the RF carrier. 
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Table 4-1. List of Figures Showing Pulse and Bar Responses 
of Specified FM Systems 



♦Note: B for Buttersworth; C for Chebyshev 
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Table 4-2. Half-Power Bandwidths Corresponding to the 
0. 1-dB Ripple Bandwidths for Three- and Four-Pole Chebyshev Filters 


Ripple 

Bandwidth 

(MHz) 

3 Poles 

4 Poles 

10 

13. 89 

12.13 

12 

16.67 

14. 55 

15 

20. 84 

18.20 

20 

— 

24.26 

30 

— 

36.39 
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Figure 4-1. T- Pulse Test Signal 
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Figure 4-2. T-Bar Test Signal 
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Figure 4-3. Response of a Flat FM Television System 
Having a 3-Pole Butterworth Filter to a T-Pulse Test Signal 
(B = 10 MHz, D p = 5 MHz) 
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Figure 4-4. Response of a Flat FM Television System 

Having a 3-Pole Butterworth Filter to a T-Bar Test Signal 

(B = 10 MHz, D = 5 MHz) 
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Figure 4-5. Response of a Flat FM Television System 
Having a 3-Pole Butterworth Filter to a T-Pulse Test Signal 
(B = 12 MHz, D - 5 MHz) 
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Figure 4-6. Response of a Flat FM Television System 
Having a 3-Pole Butterworth Filter to a T-Bar Test Signal 
(B = 12 MHz, Dp = 5 MHz) 


OUTPUT (VOLTS) 



Figure 4-7. Response of a Flat FM Television System 
Having a 3-Pole Butterworth Filter to a T-Pulse Test Signal 
(B = 15 MHz, D p = 5 MHz) 
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Figure 4-8. Response of a Flat FM Television System 

Having a 3-Pole Butterworth Filter to a T-Bar Test Signal 

(B =15 MHz, D = 5 MHz) 
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Figure 4-9. Response of a Flat FM Television System 
Having a 4-Pole Butterworth Filter to a T-Pulse Test Signal 
(B = 10 MHz, D = 5 MHz) 


4-17 


OUTPUT (VOLTS) 



TIME (H SEC) 


Figure 4-10. Response of a Flat FM Television System 
Having a 4-Pole Butterworth Filter to a T-Bar Test Signal 
(B = 10 MHz, D = 5 MHz) 
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Figure 4-11. Response of a Flat FM Television System 
Having a 4-Pole Butter worth Filter to a T-Pulse Test Signal 
(B = 12 MHz, D p = 5 MHz) 
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Figure 4-12. Response of a Flat FM Television System 
Having a 4-Pole Butterworth Filter to a T-Bar Test Signal 
(B = 12 MHz, D p = 5 MHz) 
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Figure 4-13. Response of a Flat FM Television System 

Having a 4-Pole Butterworth Filter to a T-Pulse Test Signal 

(B = 15 MHz, D = 5 MHz) 
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Figure 4-14. Response of a Flat FM Television System 
Having a 4-Pole Butterworth Filter to a T-Bar Test Signal 
(B = 15 MHz, D p = 5 MHz) 
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Figure 4-15. Response of a Flat FM Television System 
Having a 4-Pole Butter worth Filter to a T- Pulse Test Signal 
(B = 20 MHz, D p = 5 MHz) 
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Figure 4-16. Response of a Flat FM Television System 
Having a 4-Pole Butterworth Filter to a T-Bar Test Signal 
(B = 20 MHz, D = 5 MHz) 
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Figure 4-17. Response of a Flat FM Television System 
Having a 4-Pole Butterworth Filter to a T-Pulse Test Signal 
(B = 30 MHz, D p » 5 MHz) 
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Figure 4-18. Response of a Flat FM Television System 
Having a 4-Pole Butterworth Filter to a T-Bar Test Signal 
<B = 30 MHz, 0^ - 5 MHz) 
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Figure 4-19. Response of a Flat FM Television System 
Having a 3-Pole Chebyshev Filter to a T-Pulse Test Signal 

(B = 10 MHz, D = 5 MHz, Ripple = 0. 1 dB) 
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Figure 4-20. Response of a Flat FM Television System 

Having a 3-Pole Chebyshev Filter to a T-Bar Test Signal 

(B = 10 MHz, D = 5 MHz, Ripple - 0. 1 dB) 
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Figure 4-21. Response of a Flat FM Television System 
Having a 3-Pole Chebyshev Filter to a T-Pulse Test Signal 

(B = 12 MHz, D = 5 MHz, Ripple = 0. 1 dB) 
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Figure 4-22. Response of a Flat FM Television System 

Having a 3-Pole Chebyshev Filter to a T-Bar Test Signal 

(B = 12 MHz, D = 5 MHz, Ripple = 0.1 dB) 
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Figure 4-23. Response of a Flat FM Television System 
Having a 3-Pole Chebyshev Filter to a T-Pulse Test Signal 
(B = 15 MHz, D = 5 MHz, Ripple = 0. 1 dB) 
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Figure 4-24, Response of a Flat FM Television System 

Having a 3-Pole Chebyshev Filter to a T-Bar Test Signal 

(B = 15 MHz, D = 5 MHz, Ripple = 0.1 dB) 
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Figure 4-25. Response of a Flat FM Television System 
Having a 4-Pole Chebyshev Filter to a T-Puise Test Signal 
(B = 10 MHz, D = 5 MHz, Ripple = 0. 1 dB) 
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Figure 4-26. Response of a Flat FM Television System 

Having a 4-Pole Chebyshev Filter to a T-Bar Test Signal 

(B = 10 MHz, D = 5 MHz, Ripple = 0.1 dB) 
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Figure 4-27. Response of a Flat FM Television System 
Having a 4-Pole Chebyshev Filter to a T-Pulse Test Signal 
(B » 12 MHz, D = 5 MHz, Ripple -0.1 dB) 
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Figure 4-28. Response of a Flat FM Television System 
Having a 4- Pole Chebyshev Filter to a T-Bar Test Signal 
(B = 12 MHz, D = 5 MHz, Ripple = 0. 1 dB) 
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Figure 4-29. Response of a Flat FM Television System 
Having a 4-Pole Chebyshev Filter to a T-Pulse Test Signal 
(B = 15 MHz, Dp = 5 MHz, Ripple = 0.1 dB) 
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Figure 4-30. Response of a Flat FM Television System 
Having a 4-Pole Chebyshev Filter to a T-Bar Test Signal 
(B - 15 MHz, D = 5 MHz, Ripple = 0. 1 dB) 
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Figure 4-31. Response of a Flat FM Television System 
Having a 4-Pole Chebyshev Filter to a T-Pulse Test Signal 
(B = 20 MHz, D p = 5 MHz, Ripple = 0. 1 dB) 
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Figure 4-32. Response of a Flat FM Television System 
Having a 4-Pole Chebyshev Filter to a T-Bar Test Signal 
(B = 20 MHz, D = 5 MHz, Ripple = 0. 1 dB) 
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Figure 4-33. Response of a Flat FM Television System 
Having a 4-Pole Chebyshev Filter to a T-Pulse Test Signal 
(B = 30 MHz, D p = 5 MHz, Ripple -0.1 dB) 
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Figure 4-34. Response of a Flat FM Television System 
Having a 4-Pole Chebyshev Filter to a T-Bar Test Signal 
(B = 30 MHz, D p = 5 MHz, Ripple = 0. 1 dB) 
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Figure 4-35. Response of a Flat FM Television System 
Having a 3 -Pole Butte rworth Filter to a T-Pulse 
Test Signal (B = 10 MHz, D p = 6 MHz) 
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Figure 4-36. Response of a Flat FM Television System 
Having a 3-Pole Butterworth Filter to a T-Bar Test Signal 
(B = 10 MHz, D p = 6 MHz) 
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Figure 4-37. Response of a Flat FM Television System 
Having a 3-Pole Butterworth Filter to a T-Pulse 
Test Signal (B = 12 MHz, Dp = 6 MHz) 
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Figure 4-38. Response of a Flat FM Television System 
Having a 3-Pole Butterworth Filter to a T-Bar Test Signal 
(B = 12 MHz, D p = 6 MHz) 
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Figure 4-39. Response of a Flat FM Television System 

Having a 3 -Pole Butterworth Filter to a T- Pulse 

Test Signal (B = 15 MHz, D =6 MHz) 
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Figure 4-40. Response of a Flat FM Television System 
Having a 3-Pole Butterworth Filter to a T-Bar Test Signal 
(B = 15 MHz, D = 6 MHz) 
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Figure 4-41. Response of a Flat FM Television System 
Having a 4-Pole Butterworth Filter to a T-Pulse 
Test Signal (B = 10 MHz, Dp = 6 MHz) 
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Figure 4-42. Response of a Flat FM Television System 
Having a 4-Pole Butterworth Filter to a T-Bar Test Signal 
(B = 10 MHz, D = 6 MHz) 
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Figure 4-43. Response of a Flat FM Television System 
Having a 4-Pole Butterworth Filter to a T -Pulse 
Test Signal (B = 12 MHz, ~ 6 MHz) 
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Figure 4-44. Response of a Flat FM Television System 
Having a 4-Pole Butterworth Filter to a T-Bar Test Signal 
(B = 12 MHz, D p = 6 MHz) 
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Figure 4-45. Response,of a Flat FM Television System 
Having a 4-Pole Bntterworth Filter to a T-Pulse 
Test Signal (B = 15 MHz, D p - 6 MHz) 
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Figure 4-46. Response of a Flat FM Television System 
Having a 4-Pole Butterworth Filter to a T-Bar Test Signal 
(B = 15 MHz, D p = 6 MHz) 
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Figure 4-47. Response of a Flat FM Television System 
Having a 4-Pole Butterworth Filter to a T-Pulse 
Test Signal (B = 20 MHz, D - 6 MHz) 


4-55 


OUTPUT (VOLTS) 


0.7 


0.6 

0.5 

0.4 

0.3 

0.2 

0.1 

0 

- 0.1 

- 0.2 

-0.3 




0 0.1 0.2 0.3 25.0 25. 1 25.2 25.3 25.4 

TIME (fX SEC) 


Figure 4-48. Response of a Flat FM Television System 
Having a 4-Pole Butterworth Filter to a T-Bar Test Signal 
(B = 20 MHz, 0^ - 6 MHz) 
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Figure 4-49. Response of a Flat FM Television System 
Having a 4-Pole Butterworth Filter to a T-Pulse 
Test Signal (B = 30 MHz, = 6 MHz) 
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Figure 4-50. Response of a Flat FM Television System 
Having a 4-Pole Butterworth Filter to a T-Bar Test Signal 
(B = 30 MHz, D = 6 MHz) 
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Figure 4-51. Response of a Flat FM Television System Having 
a 3-Pole Chebyshev Filter to a T-Pulse Test Signal 
(B = 10 MHz, D p = 6 MHz, Ripple = 0.1 dB) 
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Figure 4-52. Response of a Flat FM Television System 
Having a 3-Pole Chebyshev Filter to a T-Bar Test Signal 
(B = 10 MHz, D p - 6 MHz, Ripple = 0. 1 dB) 
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Figure 4-53. Response of a Flat FM Television System Having 
a 3-Pole Chebyshev Filter to a T-Pulse Test Signal 
(B = 12 MHz, D = 6 MHz, Ripple = 0.1 dB) 
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Figure 4-54. Response of a Flat FM Television System 
Having a 3-Pole Chebyshev Filter to a T-Bar Test Signal 
(B = 12 MHz, D p = 6 MHz, Ripple = 0. 1 dB) 



Figure 4-55. Response of a Flat FM Television System Having 

a 3-Pole Chebyshev Filter to a T-Pulse Test Signal 

(B - 15 MHz, D = 6 MHz, Ripple = 0. 1 dB) 
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Figure 5-56. Response of a Flat FM Television System 
Having a 3-Pole Chebyshev Filter to a T-Bar Test Signal 
(B = 15 MHz, D p = 6 MHz, Ripple = 0.1 dB) 
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Figure 4-57. Response of a Flat FM Television System Having 
a 4-Pole Chebyshev Filter to a T-Pulse Test Signal 

(B = 10 MHz, D = 6 MHz, Ripple = 0. 1 dB) 
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Figure 4-58. Response of a Flat FM Television System 
Having a 4-Pole Chebyshev Filter to a T-Bar Test Signal 
(B = 10 MHz, D p = 6 MHz, Ripple = 0. 1 dB) 
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Figure 4-59. Response of a Flat FM Television System Having 
a 4-Pole Chebyshev Filter to a T-Pulse Test Signal 
(B = 12 MHz, D = 6 MHz, Ripple = 0. 1 dB) 
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Figure 4-60. Response of a Flat FM Television System 
Having a 4-Pole Chebyshev Filter to a T-Bar Test Signal 
(B = 12 MHz, D = 6 MHz, Ripple = 0. 1 dB) 
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Figure 4-61. Response of a Flat FM Television System Having 

a 4-Pole Chebyshev Filter to a T-Pulse Test Signal 

(B = 15 MHz, D - 6 MHz, Ripple = 0.1 dB) 

P 


4-70 




.0 25. 1 25.2 25.3 25.4 


2 (/i SEC) 

FM Television System 
to a T-Bar Test Signal 
Ripple = 0.1 dB) 


OUTPUT (VOLTS) 



TIME (jl SEC) 


Figure 4-63. Response of a Flat FM Television System Having 
a 4- Pole Chebyshev Filter to a T-Pulse Test Signal 
(B = 20 MHz, D = 6 MHz, Ripple = 0.1 dB) 
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Figure 4-64. Response of a Flat FM Television System 
Having a 4-Pole Chebyshev Filter to a T-Bar Test Signal 
(B = 20 MHz, D p = 6 MHz, Ripple = 0. 1 dB) 
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Figure 4-65. Response of a Flat FM Television System Having 
a 4-Pole Chebyshev Filter to a T-Pulse Test Signal 
(B =30 MHz, Dp - 6 MHz, Ripple = 0. 1 dB) 
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Figure 4-66. Response of a Flat FM Television System 
Having a 4-Pole Chebyshev Filter to a T-Bar Test Signal 
(B = 30 MHz, D = 6 MHz, Ripple = 0. 1 dB) 
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Figure 4-67. Response of a Preemphasized FM Television System 

Having a 3-Pole Butterworth Filter to a T-Pulse Test Signal 

(B = 10 MHz, D = 5 MHz) 
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Figure 4-68. Response of a Preemphasized FM Television System 

Having a 3-Pole Butterworth Filter to a T-Bar Test Signal 

(B = 10 MHz, D = 5 MHz) 
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Figure 4-69. Response of a Preemphasized FM Television System 
Having a 3-Pole Butterworth Filter to a T-Pulse Test Signal 
(B = 12 MHz, D p = 5 MHz) 
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Figure 4-70. Response of a Preemphasized FM Television System 
Having a 3-Pole Butterworth Filter to a T-Bar Test Signal 
(B = 12 MHz, D p = 5 MHz) 
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Figure 4-71. Response of a Preemphasized FM Television System 
Having a 3-Pole Butterworth Filter to a T-Pulse Test Signal 
(B = 15 MHz, Dp = 5 MHz) 
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Figure 4-72. Response of a Preemphasized FM Television System 
Having a 3 -Pole Butterworth Filter to a T-Bar Test Signal 
(B = 15 MHz, D p = 5 MHz) 
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Figure 4-73. Response of a Preemphasized FM Television System 
Having a 4-Pole Butterworth Filter to a T-Pulse Test Signal 
(B = 10 MHz, D = 5 MHz) 
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Figure 4-74. Response of a Preemphasized FM Television System 
Having a 4-Pole Butterworth Filter to a T-Bar Test Signal 
(B = 10 MHz, D p = 5 MHz) 
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Figure 4-75. Response of a Preemphasized FM Television System 
Having a 4-Pole Butterworth Filter to a T-Pulse Test Signal 
(B = 12 MHz, D = 5 MHz) 
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Figure 4-76. Response of a Preeamphasized FM Television System 
Having a 4-Pole Butterworth Filter to a T-Bar Test Signal 
(B = 12 MHz, D p = 5 MHz) 
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Figure 4-77. Response of a Preemphasized FM Television System 

Having a 4-Pole Butterworth Filter to a T-Pulse Test Signal 

(B = 15 MHz, D = 5 MHz) 
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Figure 4-79. Response of a Preemphasized FM Television System 
Having a 4- Pole Butterworth Filter to a T-Pulse Test Signal 
(B = 20 MHz, D = 5 MHz) 
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Figure 4-80. Response of a Preemphasized FM Television System 
Having a 4-Pole Butterworth Filter to a T-Bar Test Signal 
(B = 20 MHz, D = 5 MHz) 
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Figure 4-81. Response of a Preemphasized FM Television System 

Having a 4 -Pole Butterworth Filter to a T- Pulse Test Signal 

(B = 30 MHz, D = 5 MHz) 
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Figure 4-82. Response of a Preemphasized FM Television System 
Having a 4- Pole Butterworth Filter to a T-Bar Test Signal 
(B = 30 MHz, D p = 5 MHz) 
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Figure 4-83. Response of a Preemphasized FM Television System 
Having a 3-Pole Chebyshev Filter to a T-Pulse Test Signal 
(B = 10 MHz, D p - 5 MHz, Ripple = 0. 1 dB) 
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Figure 4-84. Response of a Preemphasized FM Television System 
Having a 3-Pole Chebyshev Filter to a T-Bar Test Signal 
(B = 10 MHz, D p = 5 MHz, Ripple = 0.1 dB) 
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Figure 4-85. Response of a Preemphasized FM Television System 
Having a 3-Pole Chebyshev Filter to a T-Pulse Test Signal 
(B = 12 MHz, D p = 5 MHz, Ripple = 0. 1 dB) 
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Figure 4-86. Response of a Preemphasized FM Television System 
Having a 3-Pole Chebyshev Filter to a T-Bar Test Signal 
(B = 12 MHz, D = 5 MHz, Ripple = 0. 1 dB) 
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Figure 4-87. Response of a Preemphasized FM Television System 
Having a 3 -Pole Chebyshev Filter to a T-Pulse Test Signal 
(B = 15 MHz, D p = 5 MHz, Ripple = 0. 1 dB) 
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Figure 4-88. Response of a Preemphasized FM Television System 

Having a 3-Pole Chebyshev Filter to a T-Bar Test Signal 

(B = 15 MHz, D = 5 MHz, Ripple = 0. 1 dB) 
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Figure 4-89. Response of a Preemphasized FM Television System 

Having a 4-Pole Chebyshev Filter to a T-Pulse Test Signal 

(B = 10 MHz, D = 5 MHz, Ripple = 0. 1 dB) 
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Figure 4-90. Response of a Preemphasized FM Television System 
Having a 4-Pole Chebyshev Filter to a T-Bar Test Signal 
(B = 10 MHz, D p - 5 MHz, Ripple = 0. 1 dB) 
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Figure 4-91. Response of a Preemphasized FM Television System 
Having a 4-Pole Chebyshev Filter to a T-Pulse Test Signal 
(B = 12 MHz, D = 5 MHz, Ripple = 0.1 dB) 
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Figure 4-92. Response of a Preemphasized FM Television System 

Having a 4-Pole Chebyshev Filter to a T-Bar Test Signal 

(B = 12 MHz, D =5 MHz, Ripple = 0. 1 dB) 
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Figure 4-93. Response of a Preemphasized FM Television System 
Having a 4-Pole Chebyshev Filter to a T-Pulse Test Signal 
(B = 15 MHz, D = 5 MHz, Ripple = 0.1 dB) 
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Figure 4-94. Response of a Preemphasized FM Television System 
Having a 4-Pole Chebyshev Filter to a T-Bar Test Signal 
(B = 15 MHz, D = 5 MHz, Ripple = 0. 1 dB) 
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Figure 4-95. Response of a Preemphasized FM Television System 
Having a 4-Pole Chebyshev Filter to a T-Pulse Test Signal 
(B = 20 MHz, D =5 MHz, Ripple = 0. 1 dB) 
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Figure 4-96. Response of a Preemphasized FM Television System 
Having a 4-Pole Chebyshev Filter to a T-Bar Test Signal 
(B = 20 MHz, D = 5 MHz, Ripple = 0. 1 dB) 
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Figure 4-97. Response of a Preemphasized FM Television System 
Having a 4-Pole Chebyshev Filter to a T-Pulse Test Signal 
(B = 30 MHz, D =5 MHz, Ripple = 0.1 dB) 
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Figure 4-98. Response of a Preemphasized FM Television System 
Having a 4-Pole Chebyshev Filter to a T-Bar Test Signal 
(B = 30 MHz, D p = 5 MHz, Ripple = 0. 1 dB) 
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Figure 4-99. Response of a Preemphasized FM Television System 
Having a 3-Pole Butterworth Filter to a T-Pulse Test Signal 
(B = 10 MHz, D p = 6 MHz) 
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Figure 4-100. Response of a Preemphasized FM Television System 

Having a 3-Pole Butterworth Filter to a T-Bar Test Signal 

(B = 10 MHz, D = 6 MHz) 
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Figure 4-101. Response of a preemphasized FM Television System 
Having a 3-Pole Butterworth Filter to a T-Pulse Test Signal 
(B = 12 MHz, D = 6 MHz) 
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Figure 4-102. Response of a Preemphasized FM Television System 
Having a 3 -Pole Butterworth Filter to a T-Bar Test Signal 
(B = 12 MHz, D = 6 MHz) 
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Figure 103. Response of a Preemphasized FM Television System 
Having a 3-Pole Butterworth Filter to a T-Pulse Test Signal 
(B = 15 MHz, D p * 6 MHz) 
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Figure 4-104. Response of a Preemphasized FM Television System 
Having a 3-Pole Butterworth Filter to a T-Bar Test Signal 
(B = 15 MHz, D p = 6 MHz) 
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Figure 4-105. Response of a Preemphasized FM Television System 
Having a 4-Pole Butterworth Filter to a T-Pulse Test Signal 
(B = 10 MHz, D p = 6 MHz) 
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Figure 4-106. Response of a Preemphasized FM Television System 
Having a 4- Pole Butterworth Filter to a T-Bar Test Signal 
(B = 10 MHz, D p = 6 MHz) 
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Figure 4-107. Response of a Preemphasized FM Television System 

Having a 4-Pole Butterworth Filter to a T-Pulse Test Signal 

(B = 12 MHz, D = 6 MHz) 
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Figure 4-108. Response of a Preemphasized FM Television System 

Having a 4- Pole Butterworth Filter to a T-Bar Test Signal 

(B = 12 MHz, D - 6 MHz) 
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Figure 4-109. Response of a Preemphasized FM Television System 
Having a 4-Pole Butterworth Filter to a T-Pulse Test Signal 
(B = 15 MHz, D = 6 MHz) 
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Figure 4-110. Response of a Preemphasized FM Television System 
Having a 4-Pole Butterworth Filter to a T-Bar Test Signal 
(B = 15 MHz, D p = 6 MHz) 


OUTPUT (VOLTS) 



TIME (fl SEC) 


Figure 4-111. Response of a Preemphasized FM Television System 
Having a 4-Pole Butterworth Filter to a T-Pulse Test Signal 
(B = 20 MHz, D p = 6 MHz) 
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Figure 4-112. Response of a Preemphasized FM Television System 

Having a 4-Pole Butterworth Filter to a T-Bar Test Signal 

(B = 20 MHz, D = 6 MHz) 
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Figure 4-113. Response of a Preemphasized FM Television System 

Having a 4-Pole Butterworth Filter to a T-Pulse Test Signal 

(B = 30 MHz, D = 6 MHz) 
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Figure 4-114. Response of a Preemphasized FM Television System 
Having a 4-Pole Butterworth Filter to a T-Bar Test Signal 
(B = 30 MHz, D p = 6 MHz) 
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Figure 4-115. Response of a Preemphasized FM Television System 
Having a 3-Pole Chebyshev Filter to a T-Pulse Test Signal 
(B = 10 MHz, Dp = 6 MHz, Ripple = 0. 1 dB) 
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Figure 4-116. Response of a Preemphasized FM Television System 
Having a 3-Pole Chebyshev Filter to a T-Bar Test Signal 
(B - 10 MHz, D = 6 MHz, Ripple = 0. 1 dB) 
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Figure 4-117. Response of a Preemphasized FM Television System 
Having a 3-Pole Chebyshev Filter to a T-Pulse Test Signal 

(B = 12 MHz, D = 6 MHz, Ripple = 0.1 dB) 
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Figure 4-118. Response of a Preemphasized FM Television System 

Having a 3-Pole Chebyshev Filter to a T-Bar Test Signal 

(B = 12 MHz, D = 6 MHz, Ripple = 0.1 dB) 
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Figure 4-119. Response of a Preemphasized FM Television System 

Having a 3-Pole Chebyshev Filter to a T-Pulse Test Signal 

(B = 15 MHz, D = 6 MHz, Ripple = 0. 1 dB) 
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Figure 4-120. Response of a Preemphasized FM Television System 
Having a 3-Pole Chebyshev Filter to a T-Bar Test Signal 

(B = 15 MHz, D = 6 MHz, Ripple = 0. 1 dB) 
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Figure 4-121. Response of a Preemphasized FM Television System 
Having a 4-Pole Chebyshev Filter to a T-Pulse Test Signal 
(B = 10 MHz, Dp - 6 MHz, Ripple = 0.1 dB) 
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Figure 4-122. Response of a Preemphasized FM Television System 
Having a 4-Pole Chebyshev Filter to a T-Bar Test Signal 
(B = 10 MHz, D = 6 MHz, Ripple = 0. 1 dB) 
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Figure 4-123. Response of a Preemphasized FM Television System 
Having a 4-Pole Chebyshev Filter to a T-Pulse Test Signal 
(B = 12 MHz, D p = 6 MHz, Ripple = 0. 1 dB) 
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Figure 4-124. Response of a Preemphasized FM Television System 
Having a 4-Pole Chebyshev Filter to a T-Bar Test Signal 
(B = 12 MHz, D = 6 MHz, Ripple = 0. 1 dB) 
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Figure 4-125. Response of a Preemphasized FM Television System 
Having a 4-Pole Chebyshev Filter to a T-Pulse Test Signal 
(B = 15 MHz, D p = 6 MHz, Ripple = 0. 1 dB) 
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Figure 4-126 Response of a Preemphasized FM Television System 
Having a 4-Pole Chebyshev Filter to a T-Bar Test Signal 
<B = 15 MHz, D = 6 MHz, Ripple = 0. 1 dB) 
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Figure 4-127. Response of a Preemphasized FM Television System 
Having a 4-Pole Chebyshev Filter to a T-Pulse Test Signal 
(B = 20 MHz, D = 6 MHz, Ripple = 0. 1 dB) 
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Figure 4-128. Response of a Preemphasized FM Television System 
Having a 4-Pole Chebyshev Filter to a T-Bar Test Signal 
(B = 20 MHz, D = 6 MHz, Ripple = 0. 1 dB) 












OUTPUT (VOLTS) 



Figure 4-129. Response of a Preemphasized FM Television System 
Having a 4-Pole Chebyshev Filter to a T-Pulse Test Signal 
(B * 30 MHz, D = 6 MHz, Ripple = 0. 1 dB) 
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Figure 4-130. Response of a Preemphasized FM Television System 
Having a 4- Pole Chebyshev Filter to a T-Bar Test Signal 
(B = 30 MHz, D = 6 MHz, Ripple = 0. 1 dB) 
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Figure 4-131. Response of a 525-Line Preemphasis Filter 
to a T-Pulse (T = 100 ns) 
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Figure 4-132. Response of a 525-Line Preemphasis Filter 
to the Rise Portion of a T-Bar (T = 100 ns) 
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SECTION 5 - CONCLUSIONS 


Several conclusions can be drawn from this frequency modulation distortion analysis 
for television transmission. The conclusions are based on the computed distortions 
of a T-pulse-and-bar test signal which passes through the computer-simulated system. 
The conclusions are: 

An FM television transmission system can be modeled for computer simulation of a 
well defined signal in the time domain from the video input terminals to the video 
output terminals. 

If preemphasis is employed, a Chebyshev predetection filter produces more signal 
distortion than a Butterworth filter having four poles and the same half-power band- 
width. Without preemphasis no significant difference between Butterworth and 
Chebyshev filters was observed. 

Whether a Butterworth or Chebyshev filter is used, a four-pole filter produces more 
distortion than a three-pole filter. 

An increase in the RF predetection bandwidth reduces the system distortion if all 
other parameters remain fixed. 

An increase in the peak deviation of a system increases the system distortion if all 
other parameters remain fixed. 

The T-pulse-and-bar test signal is distorted more in a system utilizing preemphasis 
than in a flat system, which does not employ preemphasis. 

The case producing the worst distortion (Figure 4-121) is a preemphasized system 
utilizing a four-pole Chebyshev predetection filter, having a ripple bandwidth of 
10.0 MHz, corresponding to a half-power bandwidth of 12. 13 MHz. This bandwidth, 
which is approximately twice the peak deviation of 6.0 MHz and significantly less 
than Carson's bandwidth, produces an overshoot distortion of only 9.9%-less than 
the C.C.I. R. limit of 13% (Reference 6). This is achieved without any equalization 
in the system. 
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APPENDIX A - PREEMPHASIS AND DEEMPHASIS FILTER ANALYSIS 


A. 1 INTRODUCTION 


The impulse responses of preemphasis and deemphasis filters for 525-line FM television 
transmission are derived from the circuits recommended by C.C. I. R. (Reference 1). 

The transfer functions are first derived from the specified circuits. Then the impulse 
responses are calculated from the transfer functions. 

In Reference 1 the configurations and component values of preemphasis and deemphasis 
filters are specified for different television systems. For the case of the 525-line tele- 
vision system, the filter configurations are shown in Figure A-l and A-2, while their 
component values are tabulated in Tables A-l and A-2, respectively. 

Figure A-3 shows the preemphasis magnitude and phase characteristics for a 525-line 
television system which, besides the intrinsic frequency response of a passive filter, 
includes a constant amplification of 3.39000 dB. This is equivalent to a voltage multi- 
plication factor of 1.479. Similarly for the deemphasis characteristic, a voltage ampli- 
fication factor of 3. 162278 is included. 

A. 2 DETERMINATION OF PREEMPHASIS AND DEEMPHASIS TRANSFER FUNCTIONS 

Consider a circuit which consists of a generator at the input terminals and a black-box 
network coupled to a load at the output terminals, as shown in Figure A-4. If the black- 
box network is a T-section, as shown in the same figure, the elements of the impedance 
matrix can be written in terms of the branch impedances of the T-section as follows: 


Z 11 = Z a + Z c ^ 
Z 12 = Z c / 

Z 22 ” Z b + Z c ) 


(A-l) 


A-l 




If V and V denote the input and output voltages respectively, if az is defined by 

1 ^ Q 

iz = z z -z , and if Z denotes the load at the output terminals, then the transfer 
11 22 12 L 

function, H is), defined as the ratio of V to V , is given by 
12 2 1 

V. 


2 Z 12 Z L 

H 12 (S) ~ V " a + z Z 
1 z 11 L 


(A-2) 


If Z , Z , Z and Z are known, the transfer function can be determined by Equation 
a b c L 

(A-2). Thus, the transfer function of both the preemphasis and deemphasis filters can 
be obtained if the filter networks are first transformed to T-sections without the parallel 
branch across the upper pair of terminals. To achieve this, the delta-star transforma- 
tion can be applied to a portion of each of the two filters. The configurations of delta and 
star networks are shown in Figure A-5 and the necessary expressions of the transforma- 


tion are given by Equation i (A-3) , 


a = 


Z 1 Z 2 


V Z 2 + Z 3 


0 = 


Z Z 
2 3 


W Z 3 


(A-3) 


y = 


Z 1 Z 3 


z + z„ + z, 

12 3 


where a, 0* and y are impedances of the arms of the star section. 


In the case of the preemphasis filter, the portion of the filter to be transformed contains 

the elements of C, R , R 0 on the left, and R on the right, which form a delta-section. 

1 z & 

Then 


1 1 + CR s 


W R 2 


(A— 4) 


A-2 



For the deemphasis filter, the portion of it involved in the transformation contains the 

elements L, R , R on the left, and R on the right. Similarly, 

1 Z Z 

LR 1 s 

Z 1 = R„ + Ls 


(A-5) 


Z„ = = R„ 

2 3 2 


From Equations (A-3) and (A -4), for the preemphasis filter, 

R 1 R 2 

R + 2R + 2R R Cs 
X ^ X ^ 

R 2 + R 1 R 2 CS 

P R + 2R + 2R R Cs 
X Z X Z 

y - a 


(A-6) 


and from Equations (A-3) and (A-5), for the deemphasis filter 


a = 

0 = 

y = 


LR iV 


2R t R 2 + (2R 2 +R 1 )Ls 


2 2 
R 1 R 2 + LR 2 S 


2R x R 2 + (2R 2 + R^Ls 


a 


(A-7) 


After a, A and y are determined, the preemphasis and deemphasis filters can be 

converted to T-sections, as shown in Figure A-6, from which the values of Z^, Z^ 

and Z and then the values of z-,„ , z .and z nn can be determined. Since Z is 75 ohms 
c 11 ±Z *Z L* 

as specified, the transfer function H (s) can be determined by Equation (A-2). 

-L te 



If H (s) and H (s) denote the transfer functions of the preemphasis and the deemphasis 
p d 


filters respectively, they can be represented by 

q -42 -11 3 

4. 57032110) + 608. 897s + 2. 12161(10) S + 2, 11953(10) s 
H ( s ) = q -4 2 -11 3 

P 2.13754(10) + 1418.22s + 3.03767(10) s + 2.11953(10) s 

1. 28 36H10) 11 + 1, 9881(10) 5 s + 8. 82881(10) 2 s 2 + 1. 02537(10) s 

H (s) = „ r 2 -8 3 

d 1.28361(10) + 2.84657(10) s + 0.205638s + 4.79565(10) s 

A. 3 DETERMINATION OF IMPULSE RESPONSES 


(A-8) 


(A-9) 


To obtain the impulse responses of the preemphasis and deemphasis filters, the inverse 


Laplace transforms of the expressions on the right side of Equations (A-8) and (A-9) 
must be obtained. These expressions are first expanded into partial fractions having the 

following form: 


H(s) = N x + N 2 


As + B 

, , 2 a 2 

(s-fij) 



(A-10) 


where N^ is a numerical constant term, 

N is the reciprocal of the coefficient of the highest-power term of the polynomial 
2 

in the denominators of the transfer functions in Equations (A-8) and (A-9), 

A, B,C are real constants to be determined, 

(c*! ± j^) is the complex root of the polynomial, 
a ^ is the real root of the polynomial. 

The roots of the polynomials are determined by a computer subroutine available in the 
CSC INFONET library. Equation (A-10) is rearranged to the form shown in Equation 
(A- 11) before taking the inverse Laplace transform. 


H(s) =N 1 + N 2 


Ats-a^ 


Aa 1 + B 


( s - a ^ + / 9 1 2 + (s-q 1 ) z + 


s-a f 


(A-ll) 


A -4 



The four terms on the right side are recognizable Laplace transforms of some elementary 
functions the transform pairs of which are listed below. 


Laplace Transforms 


Inverse Laplace Transforms 


( 1 ) 


N 


N x i(t) 


( 2 ) 


s-o. 


(s-a/ + 0* 


a l t 

e cos ( 


& 


t) u(t) 


( 3 ) 






01 l t sin (^1 


t) u(t) 


( 4 ) 


s -n 


aJ 


u{t) 


Furthermore, the second and third transforms can be combined into the form , 


K e cos ( 


0, 


t +0 ) u(t),in which K and <f> are constants to be determined. 
The partial fraction expansions of (A-8) and (A-9) are given by 

kl 


H (s) = 1 - N 
P P 


where N = 


A (s-a) 
pi 

2 k ' 0 ~Z + 


A a + 0 
P 1 P 


(s-a^+ 0 X 

1 


P 2.11953(10) 


-11 


A = 9.16115(10) 

P 


-5 


Pi 


(s-^r +0 X 


• j — + 


s-a r 


R = 503.718 

P 


C = -5.45646(10) 
P 


-9 


A-5 


(A-12) 



a = -5.4976376(10) 


± 4.3298779(10) 


a h =-3.3365331(10) 
2 


H ,(s) = 0.213813 + N, 
d d 


A d< S -“3> , A c 

1°3 

+ *d 

^3 

c a 1 

, ,2 o 2 + 

03 

(s-a 3 ) 2 + ft* s- a 4 J 


(A-13) 


where N , = 
a 


4.79565(10) 


-8 


A, = 4.43178(10) 
d 


-2 


Pd 


= 5.20966(10) 


C, = 2.31291(10) 
d 


a = -1.1754901(10) 

J 


£ = +5. 950369(10) L 


a 4 =-1.9370308(10) 

From the listed Laplace transform pairs one can easily find the time func tions for the 
terms involved in (A-12) and (A-13). Let h p (t) and h d (t) denote the respective impulse 


responses of the preemphasis and the deemphasis filters. Then 

h (t) = <5(t) - 

P 2.11953(10) 



’ V 


a. A 
2 

11 

e cos ( 


t - K 2 e 




* 


u(t) 


(A-14) 


A-6 



/? 

where = -5.4976376(10) 
1^1= 4. 3298779(10) 4 
a = -3. 3365331(10) 6 

= 9. 16262(10)"* 

<p = 1. 79422(10)” 2 

= -5.45646(10)" 9 


and 


h,(t) = 0. 213813 d(t) + 
d 


where a = -1.1754901(10) 
«3 


4.79565(10) 
6 


-8 


03 


= 5.950369(10)" 


a 4 = -1.9370308(10) 





v] 

K 3 cos ( 


t-^ 2 > + V 


u(t) 


(A-15) 


K = 0.0443185 

J 

K 4 = 2.31291(10)" 6 
4> 0 = -0. 55566(10) 

A. 4 REFERENCE 

Recommendation 405-1, Radio Relay Systems for Television (Pre-emphasis Character- 
istics for Frequency-Modulation Systems) C.C.I.R. Xllth Plenary Assembly, New Delhi, 
1970, Vol. IV, Part 1. 


A-7 






L 



Figure A-2. Deemphasis Filter for 525-Line Television 


Table A-2. Component Values of Deemphasis Filter 


L H) 

50.16 

C (p F> 

3917 

R-j mi 

275.8 

r 2 (ft) 

75 

R 3 (ft) 

20.4 
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GAIN (dB) 



PHASE (RADIANS) 








A GENERAL CIRCUIT 


T- SECTION 


Figure A-4. A General Circuit and T-Section 
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DELTA - SECTION 



Figure A-5. Delta-Section and Star-Section 


A-12 





Figure A-6. A Converted T-Section 




APPENDIX B - DERIVATION OF THE IMPULSE RESPONSE FOR 
BUTTER WORTH AND CHEBYSHEV LOW-PASS FILTERS 


B. 1 INTRODUCTION 

The impulse responses are derived from the transfer functions for two-, three-, and 
four-pole Butterworth low-pass filters and for three- and four-pole Chebyshev low-pass 
filters. The filters are assumed to have an arbitrary bandwidth of o> radians per 

-D 

second, where o> is the half-power bandwidth for the Butterworth filters and is the 
B 

ripple bandwidth for the Chebyshev filters. The general transfer functions are usually 
presented in a form such that the radian bandwidth is unity; however, replacement of 

g 

s by in those expressions results in transfer functions with a radian-frequency 

W B 

bandwidth of (o . The constants required for the Chebyshev filters are evaluated for a 
J3 

0. 1-dB ripple. 

B. 2 IMPULSE RESPONSES OF BUTTERWORTH FILTERS 


In Reference 1, the general transfer function of a low-pass Butterworth filter is given as 

h „b< s > - <B ' 1) 

n 

where B (s) is the Butterworth polynomial of order n. The coefficients of B (s), up to 
n n 

the tenth order are listed in the mentioned reference. The second, third, and fourth 
order Butterworth polynomials are given as follows: 


B (s) = s + V2 ~ s + 1 


B 3 (s) 


s 3 + 2s 2 + 2 s + 1 


(B— 2) 
(B— 3) 


B 4 (s) = s 4 + 2. 6131259s 3 + 3.4142136s 2 + 2.6131259s + 1 (B-4) 


B-l 



To obtain the transfer function of a low-pass filter with bandwidth, w D radians per 

.D 

second, the complex frequency, s, representing jw, in Equations (B-2), (B-3), and 
(B-4) is replaced by ^ . Let H 2B (s), H 3B (s), and H 4B (s) be respectively the transfer 

functions of the low-pass Butterworth filters of two, three, and four poles. Then 

2 


H 2B (S) 2 /o 2 

s + y2<w s + (o 
B B 


(B-5) 


H 3B< S > = 


3 2 2 3 

s + 2o> s" + 2(0 s 

13 13 13 


(B-6) 


H 4B< S) 


(B-7) 


s 4 + 2. 6131259&) s 3 + 3. 4142136o> 2 S 2 + 2. 6131259W 3 s + ot 
13 13 13 13 


To obtain the impulse responses, the inverse Laplace transforms of the expressions in 
Equations (B-5), (B-6), and (B-7) must be determined. The expression for H (s) in 

Lt JJ 

Equation (B-5) can be rearranged easily into a recognizable transform given by 


H 2B< S > = 


(’•I) 


(B-8) 


but the rearrangement of the other two require the use of partial fractions. The inverse 
Laplace transform of H 2B (s), given by 


h 2B (t) = 6 ^ 

is the impulse response of the two-pole filter. 


sin t u(t) 

V 2 


(B-9) 


B-2 



The expression in Equation (B-6) is expanded into partial fractions as shown in 


H 3B< S) 


A 1 (S ~ a i ] A i + B i 

+ (s 



(B-10) 


where (q^ + j^) is the complex root 


n is the real root of the denominator of (B-6) 

2 

A , B , and are coefficients of the partial fractions. 

The roots are determined by a computer subroutine in the CSC INFONET Library. 
The values of the roots and coefficients are: 

a, = -0. 

1 B 

p = + 0.8660254u> b 


ol c 


= “ U> 


A. = 


B 


B, = 0 


C =« 


B 


The corresponding impulse response , h (t), of a low-pass three-pole Butterworth filter 

o 15 

is given by 



a t a 1 



h 3B< t) = 

„ 2 „ 1 , 
C e - K e cos ( 


t " 


_ 




u(t) (B-ll) 



B-3 



Similarly for the case of the four-pole Butterworth filter, 


V 


+ B. 


H 4B (S) ' 


■ 2< y 


+ a 


where ol - -0. 92387948 w 
3 B 

3 = + 0. 38268356 ^ 

a. = -0, 38268345 w 

4 B 

/3 4 = +0. 92387951 ^ 


2 

3 



VS 

B 2 ’ 2 <“3 ' V C 2 ' D 2 


C 2 S + P 2 

2 2 2 
s - 2c* 4 s + a 4 + P 4 


(B-12) 


C 


2 



2 '“3 “V “B 


+ <“ 3 - « 4 ) ( a 3 - 3 “ 4 ) 


I ( 


- <<V' 3 4> 


+ 4(a 3 -«) 2 (a 2 +/3 2 ) 


“b ~ 2< °3~ a 4 ) °2 

<“3 '*J> ' <«4 + 4 


The two complex roots of the denominator of (B-7) are (« 3 + j/? 3 ) and (« 4 + j /3 4 ) and 
the coefficients of the partial fractions are A^, B^, C^, and D^. Since and are 
found negative, HVs) can be rearranged as 


A 2 <s 

A a + B. 
2 3 2 

|c 2 |(s-« 4 > 

| C 2 | “4 + | D 2 

4 <■-“,> * *3 

2 0 2 
< s - a 3> + S 

2 „2 
<s-a 4 ) + 



B-4 



The corresponding impulse response, h (t), is given by 

4-B 


i. 4B (t) - 


' V 

e K cos ( 

u 


t -<b 2 ) -e 


"4* 


K 3 cos ( 


‘-*3> 


u(t) 


(B-14) 




Thus, the impulse response of the Butterworth filters can be determined if the half- 
power bandwidth is specified. 

B. 3 IMPULSE RESPONSE OF CHEBYSHEV FILTERS 

The impulse responses of three- and four-pole Chebyshev filters are derived in 
this paragraph, but first some general properties of Chebyshev filters are presented. 
The frequency response of the n-pole Chebyshev filter is given by 


H nC <“> 


1 + « 2 T 2 (w) 
n 


(B-15) 


where T (w) is the Chebyshev polynomial of order n and e is the ripple factor, 
n 


B-5 



Some properties of this response are given below 


1. In the passband, the maximum value of 
T (oo) = 0, and the minimum value is 


H>) 

i nC 


is unity which occurs when 


n 


v/T 


which occurs when 


+ « 


T »<"> 


= 1 . 


2. |h ^0) | = 1 for all odd n, corresponding to the maximum values in the pass band. 

3. H JO) = ■■■■■■■■■ : for all even n, corresponding to the minimum values in the 

I nC I Vl +e 2 


pass band. 


4. 


V 1 * 


✓7 


for all n, corresponding to the minimum values in the pass band. 


+ « 


From these properties it is clear that the ratio of the maximum to the minimum power 

2 

response in the pass band is given by 1 + « . Therefore, when a ripple is specified 
by R in decibels, the ripple factor can be calculated by 

e =Vio R/1 ° - 1 (B-16) 

For example a ripple of 0.1 dB, i.e. R = 0.1, yields a ripple factor of 0.152302. 

The transfer function of a given Chebyshev filter has the form: 

G 


H „c (s > - 


n , n-1 

s + b „ s + 

n-1 


(B-17) 


+ b i s + b o 


where b . . b are positive real coefficients of the polynomial in the denominator, 
n-1 0 

and G is the real constant to be determined to satisfy the properties of H nC M 
mentioned previously. The roots, s^, of the polynomial are given in Reference 1 as 
follows: 

, (2v + 1) tt . , (2v + 1) 7T „ „ 

s = -sinh 0 sin - 1 — - — L — + j cosh <p cos — r-— 1 — ; v = 0, 1, 2, . . . , n-1 
v /n 

(B-18) 


B-6 



where sinh </> 


and cosh <f> = 


€ 2 + 1 + t j 


1 + t) 


M + 1 + T 


1/n / /l 


' 2 +1 + — 
€ e 


The constant should be determined by setting s = 0 in Equation (B-17) and noting that 


H .JO) = — — = 1, for odd n 

n<? b 0 


(B-19) 


V °> ■ r 


, for even n 


(B-20) 


as required by the properties of H Au) mentioned previously. 

For the case n = 3, the transfer function H 3C ( S ) of a low P ass three-pole Chebyshev 
filter has the form: 


H o r < s ) 9 2 2 

(s - 20^3 + + f}^) (s - a g ) 

where + j/3^ is the complex root of the polynomial, 


(B-21) 


a is the real root of the polynomial, and 

it 

is a constant to be evaluated. 

Similarly for n = 4, the transfer function H 4( _,(s) of a low-pass four-pole Chebyshev filter 
has the form: 


H 4C< S) = 


2 

(s 2 - 2 « 3 s + a 2 + /3 2 ) (s 2 - 2of 4 s + a 2 + 0 2 ) 


(B-22) 


B-7 



where (a, + j/3 ) and (a ± jy3 ) are the two complex roots of the polynomial, and 
3 3 4 4 

G is a constant to be determined. 

2 

To obtain the transfer function of a low-pass filter with ripple bandwidth , 

g 

replace s in Equations (B-21) and (B-22) by ^ g. Then 



G w 
2 B 


2 2 2 
(s^ - 2a s + £r) ( s 



(B-23) 


where a - a <o 

5 1 a 

@5 = ^1 W B 

"6 = a 2 ^ B 

and 



G w 4 
2 B 


(s 2 -2« 7 s+a^+i8^) (s 2 - 2 a g S + « g + £g) 


where a 
a 

(3 


7 

7 

8 
8 



CO 


u> 

CO 

co 


B 

B 

B 


B 


(B-24) 


To obtain the impulse responses, the right sides of Equations (B-23) and (B-24) are 
expanded into partial fractions and then the inverse Laplace transforms are taken. This 
procedure has been described previously and is not repeated here, but the results of 
calculations involved are given below. 


B-8 



The value of R is specified to be 0. 1 dB; therefore, e = 0. 152302. 

By Equation (B-18), the values of the following roots are obtained, and by Equation (B-19) 

G and G are determined: 

1 2 

« 1 = -0.484703 

P = ±1.20616 

« 2 = -0.969406 (B-25) 

G = 1. 63806 


a = -0.264156 
3 

p = ±1.12261 


«4 


= -0.637730 


8 - +0.4650 

4 

G = 0.819025 


H \ and H \ in partial 

30 K) \ B/ 

Mt) 


A i s + B i 


(s -«/ + p l 


fractions are given below. 

C. 


s - «, 


where = -C^ 


and 


B i ■ 0 


c i ■ 


- 3 

G 1 “B 


2 2 
(a -a ) ^ 

{ 5 6 ' p 5 


H 4C (.) 


A 2 S + B 2 


C 2 S + °2 




( a -« 8 ) 2 + 


(B-26) 


(B-27) 


(B-28) 


B-9 



where A = -C 
£ 


B 2 = 2(a 7 - «g) C 2 - D 2 


D, 


_ 2 (a 7 ~ “a* p 2 

2 <<?,- fs\) * <« 7 -Og) <» 7 -»“ 8) 

[(^ -j8g) + (« 7 -a 4 ) (a 7 - 3g 8 >l c 2 “b 

4 (<* 7 - a 8 ) 2 (a 2 * /3 2 ) + I (« 7 + /J,) - (“ 8 + 0j)] [(ft, ' Pg> + <“7 " 3 “g | l 


Let h (t) and h (t) denote the required impulse responses of low-pass three-pole 
3C 4C 

and four-pole Chebyshev filters with specified. Then 


h 3C (,) * 


at at 

e - e cos ( | j t - 4 > 




u(t) 


(B-29) 



(B-30) 


B-10 




Thus, the impulse responses can be determined for three- and four -pole Chebyshev 
filters if the ripple bandwidth is specified. 

B. 4 REFERENCES 


1. Louis Weinberg, Network Analysis and Synthesis , McGraw-Hill Book Company, Inc., 
1962 
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APPENDIX C - LISTING OF COMPUTER PROGRAMS 


C.l COMPUTER PROGRAM 1 


<* tpulse: to generate t-pulse 
<* dt=step size 

% 

DIMENSION X (1800) *Y < 1000) 

* 

% TIMF-INTERVALS: 

* 

T 1 = 0 . 0 
T3=0.1E-0o 
TEND=1 .8E-06 

% 

% STARTING TIME: 

% 

PT=1.0E-09 
W = 3. 14159265*5. 0E06 
T=T 1 
N= 1 

% 

500 IF ( T .GT. T3 ) GOTO 20 
A=W* ( T ) 

X ( N ) =0.70*SIN < A) **? 

Y (N) =-X (N) 

GOTO 50 
20 X ( N ) = 0 . 7 
YIN) =-0.7 
50 T = T + DT 

I F ( T .GT. TE ND ) GOTO 1000 
N=N+ 1 
GOTO 500 

* 

1000 WRITF(11**)X 
WRITE ( 12**) Y 

WRITE <6* 100) (X < I )♦ 1 = 1*1800) 
WRITE <6* 101) 

101 FORMAT (/IX* • Y ( I ) : •/) 

WRITE (6*100) (Y<I) *1=1*1800) 

100 FORMAT ( 1X*9E14 ,6 ) 

% 

STOP 

END 


C-l 



C. 2 COMPUTER PROGRAM 2 


* TPULSe: TO GENERATE T -'PULSE 
f OT=STE p SIZE 

># 

O I MENS I ON < (610) 

* “ — 

* TIME-INTERVALS: 


n = o . u 

T3=0.2E-06 
TEND=0 .HE-Ob 


* START I NO TIME: 

F 

[>T= 1 .Ok-OR 
w = 3.l4lS><326l»*i>.0F0* 
T = T 1 
M=1 


600 IF (T .GT* T 3 ) G 0 T 0 20 
A = W v T 

X (Kn =0.70*SIN< A) 

GOTO SO 
20 X(N)=0.0 
SO T = T + DT 

I F ( T .GT. TEND) GOTO 1000 
N = N* 1 
GOTO 600 


1000 WRITE ( 11 **>X 

WRITE ( 12 * 100 ) CX ( I ) « 1 = 1 ♦ 600 ) 

1 00 FORMAT ( 1 X» 9 E 14 . 6 ) 

6 

STOP 
E NO 


C. 3 COMPUTER PROGRAM 3 

+ TPULSE: TO GENERATE T-PULSE 

* DT=STEp SIZE 

% 

!' I MENS I ON X ( 19 0 0 ) 

* “ -_ 

* T IMF- INTER VALS : 

* 

T1=0. RE-06 

T2=O.7F-0b 


C-2 



T 3= 1 • 3E- 0 b 
T4 = 1 .4t -06 
TFND= 1 • 8 1 - 0 b 


^ STARTING TIME: 


GT=1 .0E-0R 

te = 3. l4lS926S#B.0E()h 

T = 0 . 0 

M=1 

400 I F ( T .GT. TDGI.TO 20 
> < M ) =-0.20 
GOTO 40 

?0 I F ( T .GT. T 2 ) G 0 f 0 30 
c = w * ( T-T l ) 

X <N) =U.7*SIN(A)**2-0.2 
GOTO 40 

30 I F ( T .GT. T 3 ) GOT 0 40 
/ (N) =0.50 
GOTO 50 

4 0 IF (T . GT. T4)(i(iro b U 
A = W* ( T-T 3 ) 

X (N)=0.70»COS<A>»»2-0.20 
GOTO 40 
GO /UM)=-().?0 
SO T = T + L> T 

JF(T .GT. Tfc.NO ) GOT 1 1 1000 
M = N + I 
GOTO SOO 


1 <1 0 0 CONTINUE; 

WRITE ( 11 .*> X 

write (6*1 on w»or*s 

1 o 1 FORMAT ( 1 X . ' T-R 4H : ERE (RAO) =*»E12.G»1X« , 0T=*»E12.6*1X»* 
• P 0 1 N T S = • * 14/) 

WRITE (6 ♦ 1 02) 

102 FORMAT ( 1 X »' OAT a STORED IN FILE.: F B A R 1 ' / / ) 

WRITE ( 10, 100) (X(I) « I = 1 « N ) 

100 FORMAT ( 1X»9E 14.6) 


STOP 

END 


C-3 



C. 4 COMPUTER PROGRAM 4 


* tpulse: to generate t -pulse 

* 0 T = S T E P SIZE 


O I MENS I ON X(l4l0) 


* TIME- INTERVALS: 


T 1 = 0 . 0 
T? = 0 • 6E-06 
Tj=0.8E-0b 
T E N D = 1 . 4E-0b 

* --------------- 

A- STARTING TIME: 

ij . 

DT = 1 .OF-CW 

•" = 3 • 14159265*5.0E08 

T = T 1 

N=1 


500 I F ( T .GT. T2)GoTO 40 
X (N) =-0 .20 
GOTO 50 

40 I F { T »6T . T3 ) GOTO 20 
A = W* < T -T<? > 

X (N) =0.70«SIN ( A) **2-0.20 
GOTO 50 
?0 X ( N ) = - 0 . 2 0 
SO T = T ♦ 0 T 

T F ( T .GT. TEND) GOTO 1000 
N = N+ 1 
GOTO 500 


1000 wRITE(ll**)X 

WRITE ( 12* 100) (X ( I ) . 1 = 1 , 1400) 

100 FORMAT UX»9£l4. 6) 


STOP 

FND 


C. 5 COMPUTER PROGRAM 5 

% PRDIMP: IMPULSE RESPONSE OF PREEMPHASIS FILTER 
DIMENSION X (2000 ) 

% 

% CONSTANTS: 

% 

A=9.161 15E-05 
B=503.718 
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C=-5.45646E-09 
ALPHA1=-S.4976376E06 
ALPHA2=-3.3365331E06 
BETA=4.3298779E04 
Al= (A*ALPHA1*B) /BETA 
A2= A 1 / A 

AK=SQRT ( A*A*A1*A1 > 

PH I = AT AN ( A2 > 

C2=1.0/2.119S3E-11 
M=2000 
T = 0 . 0 

DT=1 .0E-09 

* 

00 10 1 = 1 *N 

C1=EXP { ALPHA 1«T > *AK*COS (BETA*T-PhI ) +C*EXP ( ALPHA2°T ) 

X(I)=C1*C2 

T=T+DT 

10 CONTINUE 

% 

WRITE ( 1 1 ♦*> X 
WRITE (6*102) 

102 FORMAT (IX* • PRE I MP : PREEMPH AS I S RESPONSE '♦ 2X ♦* DAT A STORED 
IN file:fprimt»//> 

WRITE(6*101) ( X ( I ) » I = 1 ♦ N > 

101 FORMAT (1X*9E14.6) 



STOP 

END 


C. 6 COMPUTER PROGRAM 6 


* DEIMP: IMPULSE RESPONSE OF DEEMPHASIS FILTER 
DIMENSION X ( 1 80 0 ) 

% 

% CONSTANTS: 

% 

A=4.431 78F-02 
B=52096.6 
C=2«31291E-06 
ALPHAI=-1 • 1754901E06 
ALPHA2=-1 .9370308E06 
BETA=5.950369E03 
A 1 = ( A*ALPHA1+B) /BETA 
A2=A1/A 

AK = SQRT ( A*A*A1*A1 ) 

PH I = AT AN ( A2 ) 

N= 1800 
T = 0 . 0 

DT= 1 • 0E-09 
C 1 = 1 .0/4.79565E-08 
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no 10 i = 1 * n 

C2=EXP ( ALPHA 1 * T ) *AK*COS(BETA*T-PHI) + C*E XP ( ALPHA2*T ) 
X < I ) =C1 *C2 
T=T+DT 
10 CONTINUE 

% 

WRITE (6* 101 ) 

101 FORMAT <1X,*DE IMP: DEEMPHASIS RESPONSE'/) 

WRITE <6*102) 

102 FORMAT ( IX* • DATA STORED IN F I LE : FDE I MT • / ) 

WRITE <6* 103) AK*PHI *C1 

103 FORMAT (/1X»3 (E14.6*lX) /) 

WRITE < 1 1 * * ) X 

WRITE (6* 100) (X ( I ) * 1 = 1 *N> 

100 FORMAT < 1X*9E14. 6) 

% 

STOP 

END 


C. 7 COMPUTER PROGRAM 7 


* SINPUL: TO FORM COMPOSITE SIGNAL 

01 MENS I ON SIG(1800> *SIG1 (1800) *PRE<1800) *Z(10OO) *Y (1800) * * 
X (1800) ♦ W (1800) 

READ ( 1 1 * * ) S I G 
PE AD < 12 * * > S I G1 
PE AD ( 1 3 * * > PRE 

% 

Ms 1200 
DT = 1 . 0E-09 

% 

% CONVOLUTION 6. *1.479 

% 

DO 10 I = 1 * M 
S = 0.0 
51=0.0 
L = I 

no 20 j=i*i 

S=DT*SIG(L)*PRE(J)+S 
SI =DT*S I G 1 ( L ) *PRE ( J ) * S 1 
L = L-1 

20 CONTINUE 
Z < I ) = < SI G ( I >-S>*1.479 
X ( I ) = (SI61 ( I >-Sl >*1 .479 
10 CONTINUE 

» 

* FORMING THE COMPOSITE SIGNAL 

% 

4=0.2/3.162277 



AlsO, 5/3. 162277 
DO 199 I = 1 » M 
Z (I ) =Z ( I ) -A 
X ( I ) =X ( I ) ♦ A 1 

199 CONTINUE 
DO 200 1=1*600 
Y(I) = -A 

W < I ) = A 1 

200 CONTINUE 

DO 601 1=601*1800 
Y ( I ) =Z ( I-feOO) 

W ( I ) =X ( 1-600 ) 

601 CONTINUE 



100 FORMAT dX*9El4. 6) 

WR I TE (1 4 » * ) Y 
WRITE ( 15, *> W 
WRI TE (6 * 1 03) 

103 FORMAT (1 X ,» COMPOS I TE S I GN AL- 1 ♦ S TORED IN FCOMT1*/) 
WRITE(6*100) (Y(I) ,1=1*1800) 

WPITE (6* 104) 

104 FORMAT (/1X» 'COMPOSITE S I GNAL-2 » STORED IN FC0MT2 * / ) 
WRITE (6, 100) (W(I) ,1=1*1800) 

STOP 

END 


C. 8 COMPUTER PROGRAM 8 

* blNPUL! TO FORM COMPOSITE SIGNAL 
DIMENSION 816(6 01) .PRE(lHOO) , Z ( 60 1) ) ,Y (1400) 
Pb'AD( 1 1 SIG 
RE AD ( 12,«)PRE 


M=800 

PT = 1 .0E-U9 

£ 

* CONVOLUTION N *1.479 


DO 10 1 = 1 , M 

S = 0.O 

L=I 

DO 20 J=1,I 
S=DT*SIu(L> tt RrtE ( J) +S 
L =L - 1 

20 CONTINUE 
l (I) = (Sl6(I)—S)*l« i f79 
10 CONTINUE 


A = 0 » 2 / 3 * 1^)2277 
DO 199 I = 1 , M 
/ (I ) =7 (I ) -A 



199 CONTINUE 

103 FORMAT </lA. ‘COMPOSITE S I G N A L * $ T 0 R E D IN FCOMFl'//) 
DO 300 1=1*600 
Y ( I ) =-A 

300 CONTINUE 

no 301 1=601 *1400 
Y ( I ) =Z < 1-600) 

301 CONTINUE 
fc'Kl Tt ( 13**) Y 
V! KITE (6*103) 

V'PITF (6* 1 00) ( Y ( I ) * T = 1 * 1400 ) 

100 FOH-MAT (ix*9t 14. b) 

STOP 

FNO 


C.9 COMPUTER PROGRAM 9 

4 SINCOS: TO GET SIN E COS OF PHl(T) 

4 — 

4 Y ( I ) : RESULTS OF INTEGRATION 

V FLCOS: COS OF PHI(T) 

* FLSIN: SIN OF PHI (T) 

* DP: PEAK DtVl AT ION 

* 

01 MENS ION Z ( 1 40 0 ) « Y ( 1 4 0 0 ) *FLCOS(14U0) »FLSIN(1400) 
PE AO (1 1 .*) Z 

4 — 

W = 130 0 
OT=1.0E-09 
np=5. OE06 

AK=4.U*3. 14159365*pP 


•f INTEGRATION TO SET PHI(T) 


Y ( 1 ) =Z ( 1 ) *DT 
|!0 30 I =3 ♦ M 
Y ( I ) =Y ( I -1 ) +Z ( I ) *OT 
30 CONTINUE 

4 ----- — -------- — 

* TAKE COS Is SIN 


I Kl 40 I = i * M 
A A = AK* Y ( I ) 

FLCOS ( I ) =COS ( AA ) 

F LSIN ( I ) =SIN ( AA) 

40 CONTINUE 

4 

WRITE ( 12**) FLCOS 
V‘P I TE ( 13.*) FLSIN 
vRITE(6*1S0) 

iso FORMAT ( 1 X PROGRAM TO GENERATE CQS(Pril) E SlN(Phl) * 
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-- STORED IM FPC0S-.tFPSIN5'//) 

WRlTt(fe*l51)DP*DT» > 

151 FORMAT (IX. *P£ A*-OE\/I Al ION=» *El4.bf2X* 'STEP-SIZE*' ♦F12.6. * 
?X* ♦ POINTS* » ♦ 15//) 

153 FORMAT ( 1X*9E14. 6) 

WRITE <6, 155) 

155 FORMAT ( / 1X» »C;.'S( I ) : •/) 

WHITE (b» 153) (FLCOS(I) «I=1«M) 

^HITfc(6*156) 

155 FORMAT!/ 1 X . ' S I N (I ) : » // ) 

WRITE (<S* 153) (FLSIN ( I ) ♦ 1 = 1 *M) 

STOP 
F N D 


€.10 COMPUTER PROGRAM 10 


* SN3BFJ SYMMETRICAL 3-POLE BUTTERWORTH FILTER RESPONSE 

% fb: bandwidth HZ 

* DT: STEP SIZE 

* t: TIME 

* — 

DIMENSION X ( 1 8 0 0 ) 

% 

FB=15.0E06 

WB=2.0*3. 14159265*FB 
N= 1 800 
DT=1 .OE-09 
T = 0 « 0 

* 

ALPHAl=-0.5*WB 
ALPHA2=-1 .0*WB 
BETA1=0.8660254*WB 

C = WB* t " H '3/ <BETAl**2* ( ALPHA 1-ALPHA2 ) **2 ) 
B=(2*0*ALPHAl-ALPHA2)' t> C 
A = -C 

% 

AK = SQRT < A* A* ( 1 * 0 ♦ AlPH A 1 * ALPH A 1 /BET A 1 /BE T A 1 ) ) 

PHI =AT AN (ALPHA 1 /BETA 1 ) 

<* 

DO 10 I = 1 ♦ N 

X ( I ) =-AK*EXP (ALPHA 1*T) *COS ( BETA 1 *T-PHI ) ♦ CEXP ( ALPH A2*T ) 

T=T+DT 

10 CONTINUE 

0 — — -- 

WHlTE<6*100) 

100 FORMAT ( lXf • 3-POLE BUT TERWORTH * STORED IN FS3BT3 • // ) 
WRITE (6t 101 ) F B * D T ♦ n 
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101 FORMAT (IX* 

• BANDWIDTH: 

FB(HZ)='*E14.6»2X*‘STEP-SIZE=‘*% 

F 1 4 . 6 ♦ 2X 

« ‘POINTS=‘ 

*15//) 


WRITE (6* 

103) C* 

B* AK 

♦ PHI 


103 FORMAT ( IX* 

• c= • 

* E 1 4 * 6 * 

2X* 'B=* »E14.6*2X* ' AK=‘ *% 

E 1 4 .6 » 2X 

. • P H I = 

• *E1 

4.6//) 


WRITE (11 

**) X 




WR I TE ( 6 » 

102) (X 

(I) ♦ 

I = 1 » N ) 



102 FORMAT <1X*9E14.6) 

STOP 

END 


C.ll COMPUTER PROGRAM 11 

9f- SN4BF: SYMMETRICAL 4-POLE BUTTERwORTH FILTER RESPONSE 
% F6:BANDWIDTH HZ 

* OT:STEP-SIZE 

* TsTIME 

f- — 

DIMENSION X(lftOO) 

% 

FB=12.0E06 

W6=2.0*3.14159265*FB 
N=180 0 
DT=1.0E-09 
T = 0 • 0 

% 

AL PH A l=-0. 92387948 *WR 
A LPHA2=-0 .36268345* WB 
BETA1=0.38268356*WR 
6ETA2=0.92387951*WB 

% — — 

AAl=(BETAl**?-BETA2**H)+(ALPHAl-ALPHA2)*(ALPHAl-3.0*ALPHA2) 
AA2=(4.0*( ALPHA 1-AL PH A2) **2) * < ALPHA2**2 + BETA2**2) 

A A 3= ( ALPHA 1**2 + BET A 1**2) - ( ALPHA2**?+BE T A2**2 ) 

* 

0= (WB**4/(AA2+AA3*AAl) ) *AAl 

C=2 • 0* ( ALPHA 1 -ALPHA 2 ) *D/ A A 1 
B=2.0* (ALPHA 1-ALPHA2) *C-0 
A = -C 

l 

WRITE (6*100) 

100 FORMAT < IX* • 4-POLE 6UTTERWORTH*STORED IN FS4BT2*//) 

WRITE (6* 101 >FB*N*DT 

101 FORMAT (/IX* *BANDW IDTHJFB (HZ) = * » E 1 4 . 6 ♦ 2X * * PO I NTS = * ♦ % 

15* 2X* ‘STEP-SIZE (SEC) =* *E 14. 6//) 

WRITE (6* 102) ALPHA 1 , AL PHA2 ♦ BET A 1 ♦ BE T A2 

1 02 FORMAT ( /I X * • ALPHAl = ‘ *E 14 ,6*2X* • ALPHA2= • »E14.6*2A » • BETA 1 = * *% 
E14.6*2X» • BET A2= * *Kl4.6//> 

WRITE (6* 103) A*B*C*D 
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103 FORMAT (1X**A=**E14.6*2X»*B=**E14.6*2X*'C='*E14.6*2X*% 

•0=» »E14.6//> 

% 

ABC= ABS ( C > 

ABD= ABS ( 0 ) 

AK 1 =SQRT ( A*A + ( ( A* ALPH A 1 ♦ B ) /BET A 1 ) **2 ) 

PH 1 1 =AT AN ( <A*ALPHA1+B) /A/BET Al) 

AK2=SQRT ( ABC* ABC ♦ < ( ABC* ALPH A2 + ABD ) /BET A2 ) **2 ) 
PHI2=ATAN(<ABC*ALPHA2+ABD)/ABC/BETA2> 

% 

DO 10 I = 1 * N 

X(I)=AK1*EXP(ALPHA1*T)*C0S ( BET A 1 *T-PH 1 1 )-AK2*EXP(ALPHA2*T)*% 
C0S(6ETA2*T-PHI2) 

T = T + DT 
10 CONTINUE 

* “ " " 

WRITE (6*104) AK1 » PHI 1*AK2*PHI2 

1 04 FORMAT (1X*'AK1,PHI1*ETC:**4(E14.6*2X)) 

WRITE <11**>X 

WRITE (6* 105) (X ( I ) ♦ 1=1 *N> 

105 FORMAT (1X»9E14.6) 

% 

STOP 

END 


C. 12 COMPUTER PROGRAM 12 

% SN3CF: SYMMETRICAL 3-POLE CHEBYSHEV FILTER RESPONSE 
% FB! BANOWIDTH ^HZ 
% DT: STEP SIZE 

% T: TIME 

% 

DIMENSION X(IBOO) 

* 

FB = 1 5 * 0E06 

WB=2.0*3.14159265*FB 

ki = l a o o 

OT=l .OE-09 
T = 0 . 0 

4 — — 

ALPHA 1= _ 0 .4847 03*WB 
ALPHA2=-0.969406*WB 
BETAl=1.20616*WB 

C = 1 ,63806*WB**3/ (BET A 1**2+ ( ALPHA 1 -ALPH A2 ) **2) 

B= (2.0* ALPHA 1- ALPH A2) *C 
A =-C 

4 — -- 

AK = SQRT ( A* A* ( 1 .0 + ALPHA 1*ALPHA1 /BETA 1 /BETA 1 ) ) 
PHI=ATAN(ALPHA1/BETA1 ) 
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4 

DO 10 I = 1 * N 

X ( I ) =-AK*£XP (ALPHA 1*T ) *COS < BE T A 1 * T-PH I ) + C' H 'EXP(ALPHA2#T) 

T=T+DT 

10 CONTINUE 

% “ 

WRITE (6*100) 

100 FORMAT (IX* '3-POLE CHEBYSHEV * STORED IN FS3CT3 ’ // ) 

WRITE (6.101 )FB.DT»N 

101 FORMAT ( IX. • BANOW IDTHJFB ( HZ >='*E14.6*2X.'STEP-SIZE='»% 
F.14.6»2Xt *POINTS=» .15//) 

WRITE (fe»103)C»B.AK. PHI 

1 03 FORMAT (1X.»C=*.E14.6.2X»»B=»»E14.6»2X.«AK=»»* 

F 14.6* 2X» •PHI-* »E 14.6//) 

WRITE ( 1 1 .*) X 

WRITE (6.102) (X(I)*I=l*N> 

102 FORMAT ( 1X.9E14.6) 

STOP 

END 


C. 13 COMPUTER PROGRAM 13 

% SN4CF : SYMMETRICAL 4-POLE CHEBYSHEV FILTER RESPONSE 
% F8*.BANDWIDTH HZ 

* DT : STEP-S I ZE 
» T : T I ME 

% 

DIMENSION X ( 1 8 0 0 ) 

% 

FH = 30 * OE06 

WB=2.0*3. 14159265*FB 
N= 1800 
DT= 1 • OE-09 
T = 0 ♦ 0 



ALPHA1=-0.637730**R 
ALPHA2=-0.264156*WH 
BETA1=0. 4650* WB 
BETA2=1 . 12261»WB 

4 

AA1= ( BETA 1**2-BET A2**2) ♦<ALPHAl-ALPHA2)M ALPHA 1 -3 . 0* ALPHA2 ) 
AA2= (4.0* (ALPHA1-ALPHA2) **2) * < ALPH A2**2+BE T A2**2 > 

A A 3= (ALPHAl**2*BETAl-»*2>- < ALPH A2**2* BE T A2**2 ) 

* 

D= ( 0 .81 9025*1X6**4/ ( AA2+AA3*AA1 ) ) *AAl 

C=2.0* (ALPHA 1-ALPHA2) *D/AAl 
P=2 • 0* ( ALPHA 1-ALPHA2) *C-D 
A = -C 

* 

WRITE (6» 100) 
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100 FORMAT (1X**4-P0LE CHEBYSHE V * STORED IN FS4CT5 • // ) 
WRITE(6*101)FB«N*DT 

101 FORMAT (/IX* ' BANDWIDTH : FB (HZ)=»*E14.6*2X»»P0INTS=**% 

1 5 * 2X » ‘STEP-SIZE ( SEC ) = • »E14.6//) 

WRITE <6»102) ALPHA 1 . ALPHA2 ♦ BETA 1 » BETA2 

102 FORMAT (/1X» » ALPHAl=» » E 1 4 . 6 » 2X * • ALPHA2= • *E14.6*2X**B£TAl='*% 
E14.6*2X* »BETA2=» *El4.6//> 

WRITE(6*103)A*B»C*P 

103 FORMAT <1X»*A=**E14.6*2X*»B='*E14.6*2X*»C=**E14.6*2X*% 

»D=» .E14.6//) 

% — — 

ABC= ABS ( C ) 

ABD=ABS(D) 

AK1 = SQRT <A*A*((A*ALPHAl*B) /BETAl ) **2> 
PHI1=ATAN((A*ALPHA1*P)/A/BETA1) 

AK2 = SQRT ( ABC* ABC ♦ ( ( ARC*ALPHA2+ ABO ) /BET A2 ) **2 ) 

PH 1 2 = AT AN ( ( ABC* ALPHA2+ABD) /ABC/BET A2) 

% 

DO 10 I = 1 » N 

X < I ) =AK l*EXP ( ALPHA 1*T) *COS <BETAl*T-PHI 1 ) -AK2*EXP < ALPHA2*T ) ** 

COS <BETA2*T-PHI2> 

T = T* DT 
10 CONTINUE 

% 

WRITE (6* 1 04 ) AKl »PHI 1 * AK2»PHI2 

104 FORMAT < 1 X * • AK I ♦ PH 1 1 *ETC 5 » 1 4 ( E 1 4 . 6 * 2X > / ) 

WRITE (11**) X 

WRITE (6*105) (X ( I > * 1 = 1 *N) 

105 FORMAT (1X.9E14. 6) 

% 

STOP 

END 


C. 14 COMPUTER PROGRAM 14 

% CONDPH: TO GET PHID(T) 

* FLCOSU ) : COS OF PHI C T > 

* FLSIN(I): SIN OF PHI(T) 

* HBF(I): PREDETECTION FILTER RESPONSE 
% All) : A ( T ) 

% B < I ) :b (T) 

* PHI < I) :PHI <T> 

% PHID(I): DERIVATIVE OF PHI ( T ) 

% AM J MOD- 1 NDE X 
% FM:MAX MOD-FREQUENCY 

% 

DIMENSION FLCOS < 1 800 ) *FLSIN ( 1 800 ) «HBF ( 1000 ) * A < 1800 ) ♦ % 
B < 120 0) * PH 1(1200)* PH 10(1200) 

READ(ll**)FLCOS 
READ ( 12**) FLSIN 
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READ ( 1 3 ♦ * ) HBF 

% 

% constants: 

% 

M=1200 
DT= 1 • 0E-09 
DP=5.0E06 

AK=4.0*3.14159265*nP 

% - — — — 

% CONVOLUTION TO GET A(T) fc. B(T) 

% 

no 10 i = l ♦ m 

51=0.0 
S2= 0 • 0 
L = I 

DO 20 J= 1 * I 

S1=DT*HBF ( J) *FLCOS (L) +S1 
S2=DT*H8F (J)*FLSIN(L)*S2 
L = L-1 

20 CONTINUE 
A ( I ) =S 1 
P < I > =S2 
10 CONTINUE 

% 

* TO COMPUTE AT AN ( B/A) 

* 

DO 30 I = 1 * M 

I F ( A ( I ) .FQ. 0.0) GOTO 40 
C2=B ( I ) / A ( I ) 

PHI ( I ) = AT AN ( C2 > 

GOTO 30 
40 PHI ( I) =0.0 
WRITE (6*101 ) I 

101 FORMAT ( IX* ♦ A=0 .0 • *2X* • I=» ♦ 15/) 

30 CONTINUE 



% DERIVATIVE OF PKl (T) : 

% 

Nl=M-2 

DO 50 1 = 1 * N 1 

PH ID (I ) = <-3.0*PHI ( I ) ♦4.0*PHI ( I*1 ) -PHI ( I + 2) ) / (2.0*DT*AK) 

50 CONTINUE 

% 

WRITE(6*151) DP * M ♦ DT 

151 FORMAT (1X*»PEAK-DEVIATI0N: • .E14.6.2X* »P0INTS='*I5*2X,% 

* STEP-SI ZE= ' *E 1 4 . 6// ) 

VJRITE(14»*)PHID 

WRITE (6* 153) (PH ID ( I ) *1=1 *M) 

153 FORMAT ( 1X*9E14. 6) 

<* — 

STOP , 

FND 


014 



C. 15 COMPUTER PROGRAM 15 

% CONDPH: TO GET PhID(T) 

■v. FLCOS ( I ) : COS Of- PHI ( T ) 

* FLSIN ( I ) : SIN Of PHI ( T ) 

* HBF ( I ) : PREDE TEC T I ON FILTER RESPONSE 

* A ( I ) : M T ) 

* h ( I ) : h ( T ) 

* PhKDSPHUT) 

V PH I D ( I ) : DERIVATIVE OF PH I ( T ) 

* AMSMOO-INDEX 

'* fm:max MOD-FREOcEMCY 
^ — — — — — — — — — — — — — — — — — — — — — — — — — — 

DIMENSION FLCOSUaOU) «KLSIN(1400) *H8F U800) *A(1400) ♦ % 

M 1400) MHI ( 1400) MHIIJ ( 1400) ,hUE (1800)tZ<1400)*Y(l4GO) 

PEADCllf*)FLCOS 

READ ( 129*) FLSIN 

PFAD <13 9*-) HBF 

HEAD (14#*) HDE 

* -------- — 

* CONSTANTS: 

* 

M= 1200 
PT=1 .0E-09 
DP=S.UE06 

A K =4.0*3. 14159265* f'P 


* CONVOLUTION TO GET A ( T ) E. B ( T ) 

*• 

DO 10 I=i»M 
S1=0.0 
S2 = 0 • U 
L = I 

00 20 J=1»I 

Sl=DT*hBF < J)*FLCOS(L) ♦ S 1 
S2=DT*hrtF < J) *FLSIM(L) + S2 
L =L” I 

20 CONTINUE 
A < I ) =S1 
MI) =S2 
10 CONTINUE 


■y TO COMPUTE A T AN ( H / A ) 


DO 30 I=i»M 

IF ( A ( I ) . FO . U . 0 ) G ( > T 0 40 

C2 = B ( I ) /A < I ) 

PHI ( I ) = AT AN ( C2 ) 

GOTO 30 
40 PHI ( I ) =0.0 
WRITE (69 101 ) I 


C-15 



1 01 FORMAT < IX. * A=0.0 ' *2X* • 1=' * I5/> 

30 CONTINUE 

^ — — — — — — — — — — — — — — — — — — — — — — — — — 

* DERIVATIVE of PHl(T): 

> — - “ 

Nl=M-2 
no bO I=l»Nl 

phID ( I ) = (-3.0*PhI ( I ) +4.0*PHI ( I +1 ) -PHI ( 1+2) ) / (2.0*DT*AK ) 

6 0 CONTINUE 


WRITE t 6 * 1S1 ) UP.M.D7 

151 FORMAT ( IX. » Pt ak-DE V 1 AT I ON : 

♦STEP-S1ZE=» .E14.6//) 

153 FORMAT ( 1 X ♦ RE 1 A- . 6 ) 


.L14.6.2X. «P01NTS=« .I5.2X.* 


•* CORRECTION OF POINTS 

Sf — 

WRITE <6. 1 12) 

112 FORMAT ( IX. * CORRECTED PHID'/) 

1 = 1 

fiO A B = A B '3 (PhID ( I ) ) 

I F ( AB .61. 10.0)GOTO 72 
GOTO 71 

72 WRITE (6.75) I.PHID(I-I) * PH I U ( I + 2 ) 

7s FORMAT (/IX* 'ERROR PO I NT : ' ♦ 1 5 » 2X » 2E 1 4 . 6/ ) 
Al = (Ph ID ( 1-1 ) -PHIL) ( 1+2) ) /3.0 
Pri 1 0 ( I ) =PhID l 1-1 ) -A 1 
PM ID ( I + l ) =PHID ( 1-1 ) -2.0*A1 
1 = 1 + 1 

71 CONTINUE 
1 = 1 + 1 

IF (I .GT. M ) GOTO «1 
00 TO 80 

Hi WRITE HS< 153) (PHID ( I )♦ 1 = 1 *M) 

L M ^ ^ 

— — — — — — — — — — — — — — — “ 

* CONVOLUTION OF CORRECTED PHID & HDE 

* — --- 

Cl=-0. 0532562*3.162278 
N2=M-60U 

(JL 

no 201 1=601 ♦ M 

2(1-600) =PHID ( I) +0.06325 

201 CONTINUE 

DO 60 1=1 *N2 

S = 0 . 0 

L=I 

no 70 j=i«i 

S = DT*rlDE ( J) *Z (L) +s 
L=L-1 

70 CONTINUE 

Y(I)=(Z(I)*0.213813+S)*3.162278+C1 
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*0 CONTINUE 

* — — — — * 

♦KITE ( lb**) Y 
«HITE ( <S * 1 02 ) 

102 F OHMAT ( / IX, 'FINAL OUTPUT'//) 
WHITE (6* lb 3) (Y ( I ) « 1 = 1 *N2) 

STOP 

END 


C. 16 COMPUTER PROGRAM 16 

* CONDPm: TO GET PHIO(T) 

* FLCOS { I ) : COS OF PHI (T) 

S FLSIN ( I ) : SIN OF PHI (T > 

* HHF(I): pkedetection filter RESPONSE 
v a ( I ) : A ( T ) 

* B<I) SH(T) 

* PHI ( I ) tPHl ( T ) 

9 - PHI D ( I ) : DERIVATIVE OF PHl(T) 

* AM: MOO- INDEX 

<* FM1MAX MOD-FREO' iLNCY 

DIMENSION FLCOS ( 1*00 ) *FIS1 N ( lwOO ) *HBF ( 1800 )* A ( 1800 ) ♦ % 

H ( 1300) * PH I (IriOO) * hh ID ( 1800 ) » HDE (18U0)*Z(1800),Y(1600> 

READ (11**) FLCOS 

kEADI 12**>FLS1N 

PE AD ( 13**) HBF 

READ ( 14,*) HDE 


constants: 

A — — — — — — — — — — 

M = 1 2 .1 0 
DT=1 .OE-OQ 
|)P = 6.0E06 

AK = A.0*3. 1A159265*'.'P 

* CONVOLUTION TO GET A(T) $> B(T) 

Y "" 

00 10 1=1 *M 
S1=0.0 
S2=0 . 0 
1=1 

no 20 j=i*i 

S1=DT*HBF ( J) *FLCOS ( L ) ♦ S 1 
S2=DT*HBF ( J)*FLSIN(L) + S2 
L = L- 1 

20 CONTINUE 
A ( I ) =S 1 
P ( I ) =S2 
10 CONTINUE 
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+ TO COMPOTE A T A N ( B / A ) 

•F — — — 

00 3u 1 = 1 

IF (A(I) .EQ. 0.0) GOTO 40 
C2=8 ( I ) / A ( I ) 

PMI ( I ) =ATAN(C2) 

GOTO 30 
4 0 PHI ( I ) =0 . 0 
WHITE <fe» 1 a 1 ) I 

101 FORMAT (IX* ' A = 0 . 0 ' » 2 X * * 1 = ' *15/) 
30 CONTINUE 


4, DERIVATIVE OF PHI (T) : 


Ml = M-2 

DO 130 1 = 1 » Ml 

PH ID ( I ) - ( -3. 0*PriI ( I ) +4. 0*PHI ( I ♦ 1 ) -PHI ( 1+2 ) ) / (2. 0*OT*AK ) 

50 CONTINUE 

sf. — — — ----- 

WRITE (6, ibl ) DP*M,m 

151 FORMAT (IX** PEAK “DEVIATION! • * E 1 4 . 6 ♦ 2X » • PO I NTS= * * I 5 ♦ 2X » * 
•5TEP-SIZE=' ♦E14.5//) 

153 FORMAT (1X.9E14.6) 


4 CORRECTION OF POINTS 

4 --- ------ 

WRITE (60 12) 

112 FORMAT ( IX « 'CORRECTED PHID*/) 

1 = 1 

80 A8 = AS5 ( PH I D ( I ) ) 

I F ( AB .GT. 10.0) GOTO 72 
GOTO 71 

72 WRITE (6*75) I »PHID ( 1-1 ) «PHID ( 1+2) 

7b FORMAT (/ 1 X ♦» ERROR PO I NT ! * * 1 5 * 2X * 2E 1 4 . 6/ ) 
A 1 = ( PH I U ( I “ 1 > “PH 1 0 ( I ♦ 2 ) ) / 3 . 0 
PHID ( I ) =PH I D ( I — 1 ) — A 1 
PHIDU*1)=PHID(I-1)-2.0*A1 
1 = 1 + 1 

71 CONTINUE 
1 = 1 + 1 

IF ( I .GT. M) GOTO Hi 
GO TO BO 

81 wRITE(6»lb3) (PHID(I) »I=1»M) 

PP — «• — — —— — — — — — — — — — — — — ^ — — — — — — — - 

4 CONVOLUTION OF CORRECTED PHID c. HUE 

4 * ~ 

C2 = 0 .5/3. 1 522 7 8 
C 1 = 1) . bO 
N2=M-bOO 
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00 201 I=601»m 
l ( 1-600) =PHID < I ) -c? 

201 CONTINUE 
00 60 I = 1 « N2 
S= 0 • 0 
L = I 

no 70 j=i»i 
S=DT*HOE ( J) (L) +S 

L = L-1 

70 CONTINUE 

Y ( I) = (Z ( I ) *0.21381 3 + S) *3.162278 + C1 
6 0 CONTINUE 

* 

WPITE ( IS**) Y 
v. K I T E ( 6 * 10 2) 

102 FOPM<\T{/ IX* 'FINAL OUTPUT'//) 

I TE { 6 t 1 S3 ) ( Y ( I ) » I = 1 * N2 ) 

STOP 
E NO 
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APPENDIX D - GLOSSARY OF VARIABLES 
USED IN THE COMPUTER PROGRAMS 


D. 1 VARIABLES USED IN COMPUTER PROGRAM 1 

A A variable which is the product of angular frequency W and time T. 

DT An increment of time T. 

I An integer index for a DO-Loop. 

N An integer variable to count the number of points into which the time interval is 

divided. It also acts as the index of two vectors. 

T The current value of time. 

T1 A specified time interval. 

T3 A specified time interval. 

TEND The largest value of time for which calculations are made. 

W Angular frequency of a sinusoidal function. 

X A vector used to retain a set of calculated values. 

Y A vector used to retain a set of calculated values. 

D.2 VARIABLES USED IN COMPUTER PROGRAM 2 

A A variable which is the product of angular frequency W and time T. 

DT An increment of time T. 

I An integer index for a DO-Loop. 

N An integer variable to count th number of points into which the time interval is 

divided. It also acts as the index of two vectors. 

T The current value of time. 

T1 A specified time interval. 

T3 A specified time interval. 

TEND The largest value of time for which calculations are made. 

W Angular frequency of a sinusoidal function. 

X A vector used to retain a set of calculated values. 
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D. 3 VARIABLES USED IN COMPUTER PROGRAM 3 

A A variable which is the product of angular frequency W and time T. 

DT An increment of time T. 

I An integer index for a DO- Loop. 

N An integer variable to count the number of points into which the time interval is 

divided. It also acts as the index of two vectors. 

T The current value of time. 

T1 A specified time interval. 

T2 A specified time interval. 

T3 A specified time interval. 

T4 A specified time interval. 

TEND The largest value of time for which calculations are made. 

W Angular frequency of a sinusoidal function. 

X A vector used to retain a set of calculated values. 

D. 4 VARIABLES USED IN COMPUTER PROGRAM 4 

A A variable which is the product of angular frequency W and time T, 

DT An increment of time T. 

I An integer index for a DO- Loop. 

N An integer variable to count the number of points into which the time interval is 

divided. It also acts as the index of two vectors. 

T The current value of time. 

T1 A specified time interval. 

T2 A specified time interval. 

T3 A specified time interval. 

TEND The largest value of time for which calculations are made. 

W Angular frequency of a sinusoidal function. 

X A vector used to retain a set of calculated values. 
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D.5 VARIABLES USED IN COMPUTER PROGRAM 5 

A A coefficient in the partial fractions of the transfer function of the preemphasis 

filter. 

AK One of the two real constants of the impulse response of the preemphasis filter. 

ALPHA1 The real part of the complex root of the polynomial associated with the transfer 
function of the preemphasis filter. 

ALPHA2 The real root of the polynomial associated with the transfer function of the pre- 
emphasis filter. 

A1 The value of an expression used repeatedly. 

A2 The value of an expression used repeatedly. 

B A coefficient in the partial fractions of the transfer function of the preemphasis 

filter. 

BETA The imaginary part of the complex root of the polynomial associated with the 
transfer function of the preemphasis filter. 

C A coefficient in the partial fractions of the transfer function of the preemphasis 

filter . 

Ci A portion of the expression for the impulse response function of the preemphasis 

filter. 

C2 The second real constant in the impulse response function of the preemphasis 

filter. 

DT An increment of time T. 

I An integer index used in the DO- Loops. 

N An integer specifying the number of calculations. 

PHI The phase angle in the argument of the cosine term of the impulse response of 

the preemphasis filter. 

T The current value of time. 

X A vector to retain the values of the impulse response of the preemphasis filter. 
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D. 6 VARIABLES USED IN COMPUTER PROGRAM 6 

A A coefficient of the partial fractions of the transfer function of the deemphasis 

filter. 

AK One in the two real constants of the impulse response functions of the de- 

emphasis filter. 

ALPHA 1 The real part of the complex root of the polynomial associated with the transfer 
function of the deemphasis filter, 

ALPHA2 The real root of the polynomial associated with the transfer function of the de- 
emphasis filter. 

A1 The value of an expression used repeatedly. 

A 2 The value of an expression used repeatedly. 

B A coefficient of the partial fractions of the transfer function of the deemphasis 

filter. 

BETA The imaginary part of the complex root of the polynomial associated with the 
transfer function of the deemphasis filter. 

C A coefficient in the partial fractions of the transfer function of the deemphasis 

filter. 

Cl The second real constant in the impulse response function of the deemphasis 

filter. 

C2 The value of a portion of the expression for the impulse response function of 

the deemphasis filter. 

DT An increment of time T. 

I An integer index used in DO-Loops. 

N An integer specifying the number of calculations. 

PHI The phase angle in the argument of the cosine term of the impulse response of 

the deemphasis filter. 

T The current value of time. 

X A vector used to retain the sample values of the impulse response function. 
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D. 7 VARIABLES USED IN COMPUTER PROGRAM 7 

A A real constant which is used repeatedly. 

A1 A real constant used repeatedly. 

DT The value of the increment of time. 

I An index for DO- Loops. 

J An index for DO- Loops. 

L An index used to furnish a special order of vector components. 

M An integer to specify the number of calculations. 

PRE The impulse response of the preemphasis filter. 

S A variable used to form a summation in a convolution operation. 

Si A variable used to form a summation in convolution operation. 

SIG A vector to retain a test signal. 

SIG1 A vector to retain a test signal. 

X A vector to retain the results of a convolution. 

W A vector to retain the final result formed by adding two component parts of a 

test signal by the principle of superposition. 

Y A vector to retain the final result formed by adding two component parts of a 

test signal by the principle of superposition. 

Z A vector to retain the result of a convolution. 

D. 8 VARIABLES USED IN COMPUTER PROGRAM 8 

A A real constant which is used repeatedly. 

DT The value of the increment of time. 

I An index for DO -Loops. 

J An index for DO-Loops. 

L An index used to furnish a special order of vector components. 

M An integer to specify the number of calculations. 

PRE The impulse response of the preemphasis filter. 

S A variable used to form a summation in a convolution operation. 


D-5 



SIG A vector to retain a test signal. 

X A vector to retain the results of a convolution. 

Y A vector to retain the final result formed by adding two component parts 

of a test signal by the principle of superposition. 

Z A vector to retain the result of a convolution. 


D.9 VARIABLES USED IN COMPUTER PROGRAM 9 


AA 

AK 

DP 

DT 

FLCOS 

FLSIN 

I 

M 

Y 

Z 


A product of two variables used repeatedly. 

The produce of 477 times peak deviation. 

Peak frequency deviation. 

Step size of time interval used in integration. 

A vector used to retain values of the function cos 0. 

A vector used to retain values of the function sin 
An index used in DO- Loops. 

An integer specifying the number of calculations. 

A vector used to retain results of integration. 

A vector used to retain sample values of the modulating wave. 


D.10 VARIABLES USED IN COMPUTER PROGRAM 10 

A A coefficient of the partial fractions of the transfer function of a symmetrical 

3 -pole Butte rworth filter. 

AK A constant in the impulse response function of the symmetrical 3-pole Butter- 

worth filter- 

ALPHAl The real part of the complex root of the Butterworth polynomial associated 
with the transfer function of the filter. 

ALPHA2 The real root of the Butterworth polynomial associated with the transfer 
function of the filter* 

B A coefficient of the partial fractions of the transfer function of a symmetrical 

3 -pole Butterworth filter. 

BETA1 The imaginary part of the complex root of the Butterworth polynomial asso- 
ciated with the transfer function of the filter. 


DT 

FB 


Increment of time T. 

Frequency bandwidth of the filter in Hz. 
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I An integer index used in DO- Loops. 

N An integer specifying the number of calculations. 

PHI The phase angle in the argument of the cosine term of the impulse response 

function of the filter. 

T The current value of time. 

WB Angular frequency bandwidth of the filter. 

X A vector to retain sample value of the impulse response. 


D.H VARIABLES USED IN COMPUTER PROGRAM 11 


A 

AA1 

AA2 

AA3 

ABC 

ABD 

AK1 

AK2 

ALPHA1 

ALPHA2 

B 

BETA1 

BETA2 

C 

D 


A coefficient of the partial fractions of the transfer function of a symmetrical 
four-pole Butterworth filter. 

A real variable denoting an expression used repeatedly. 

A real variable denoting an expression used repeatedly. 

A real variable denoting an expression used repeatedly. 

The absolute value of C. 

The absolute value of D. 

A constant in the impulse response function of the four-pole Butterworth filter. 

A constant in the impulse response function of the four-pole Butterworth filter. 

The real part of the first complex root of the Butterworth polynomial of the 
four-pole filter. 

The real part of the second complex root of the Butterworth polynomial of the 
four -pole filter. 

A coefficient in the partial fractions of the transfer function of the four-pole 
Butterworth filter. 

The imaginary part of the first root of the Butterworth polynomial for the 
four-pole filter. 

The imaginary part of the second root of the Butterworth polynomial for the 
four-pole filter. 

A coefficient in the partial fractions of the transfer function of the four-pole 
filter. 

A coefficient in the partial fractions of the transfer function of the four-pole 
filter. 
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D. 12 VARIABLES USED IN COMPUTER PROGRAM 12 

A A coefficient in the partial fractions of the transfer function of a symmetrical 

three-pole Chebyshev filter. 

AK A constant in the impulse response function of the three-pole Chebyshev filter. 

ALPHA1 The real part of the complex root of the Chebyshev polynomial associated with 
the transfer function of the three-pole filter. 

ALPHA2 The real root of the Chebyshev polynomial associated with the transfer function 
of the three -pole filter. 

B A coefficient in the partial fractions of the transfer function of the filter. 

BETA1 The imaginary part of the complex root of the Chebyshev polynomial associated 
with the transfer function of the three-pole filter. 

C A coefficient of the partial fractions of the transfer function of the three-pole 

Chebyshev filter. 

DT Increment of time T. 

FB The frequency bandwidth of the filter. 

I An integer index used in DO-Loops. 

N An integer specifying the number of calculations. 

PHI The phase angle of the argument of the cosine term of the impulse response of 

the three-pole Chebyshev filter. 

T The current value of time. 

WB The angular frequency bandwidth of the three-pole Chebyshev filter. 

X A vector to retain the values of the impulse response of the three-pole 

Chebyshev filter. 

D. 13 VARIABLES USED IN COMPUTER PROGRAM 13 

A A coefficient in the partial fractions of the transfer function of a symmetrical 

four-pole Chebyshev filter. 

AA1 A variable representing an expression used repeatedly. 

AA2 A variable representing an expression used repeatedly. 

A A3 A variable representing an expression used repeatedly. 

ABC The absolute value of C. 

ABD The absolute value of D. 
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AK1 


AK2 

ALPHA1 

ALPHA2 

B 

BETA1 

BETA2 

C 

D 

DT 

FB 

I 

N 

PHU 

PHI2 

T 

WB 

X 


One of the two constants in the impulse response function of the four-pole 
Chebyshev filter. 

One of the two constants in the impulse response function of the four-pole 
Chebyshev filter. 

The real part of the first complex root of the Chebyshev polynomial associated 
with the transfer function of the four-pole filter. 

The real part of the second complex root of the Chebyshev polynomial associated 
with the transfer function of the four-pole filter. 

A coefficient in the partial fractions of the transfer function of the four-pole 
Chebyshev filter. 

The imaginary part of the first complex root of the Chebyshev polynomial 
associated with the transfer function of the four-pole filter. 

The imaginary part of the second complex root of the Chebyshev polynomial 
associated with the transfer function of the four-pole filter. 

A coefficient in the partial fractions of the transfer function of the four-pole 
Chebyshev filter. 

A coefficient in the partial fractions of the transfer function of the four-pole 
Chebyshev filter. 

An increment of time T. 

The frequency bandwidth of the four-pole Chebyshev filter. 

An integer index used in DO- Loops. 

An integer specifying the number of calculations. 

The phase angle in the argument of the first cosine term of the impulse response 
function of the four-pole Chebyshev filter. 

The phase angle of the second cosine term of the impulse response function of the 
the four-pole Chebyshev filter. 

The current value of time. 

The angular frequency bandwidth of the four-pole Chebyshev filter. 

A vector to retain the values of the impulse response of the four-pole Cheby- 
shev filter. 
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D. 14 VARIABLES USED IN COMPUTER PROGRAM 14 


A 

AK 

B 

C2 

DP 

DT 

FLCOS 

FLSIN 

HBF 

I 

J 

L 

M 

N1 

PHI 

PHID 

51 

52 


A vector representing the convolution of cos<£(t) with H^{t), the impulse 
response of the equivalent lowpass filter of a predetection filter. 

The constant relating peak frequency deviation to the voltage of the modulating 
wave. 

A vector representing the convolution of sin <f> (t) with hj{t). 

The quotient of a component of the vector B and a corresponding component of 
the vector A. 

The value of a specified peak frequency deviation. 

An increment of time T. 

A vector to retain sample values of cos(£(t). 

A vector to retain sample values of sin$(t). 

A vector to retain sample values of the impulse response of a prediction filter. 
An integer index used in DO- Loops. 

An integer index used in DO-Loops. 

An integer used to furnish a special ordering for vector components. 

An integer specifying the number of calculations. 

An integer specifying the number of calculations. 

A vector, the components of which are generated by the arctangent of C2. 

A vector of sample values of the derivative of PHI. 

The summation of terms involved in a convolution operation. 

The summation of terms involved in a convolution operation. 


D.I5 VARIABLES USED IN COMPUTER PROGRAM 15 

A A vector representing the convolution of cos<£»(t) with h^(t), the impulse response 

of the equivalent lowpass filter of a predetection filter. 

AB The absolute value of the derivative of PHI. 

AK The constant relating peak frequency deviation to the voltage of the modulating 

wave. 

A1 The value of an expression used repeatedly. 

B A vector representing the convolution of sin$(t) with h^t). 
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Cl 

C2 

DP 

DT 

FLCOS 

FLSIN 

HBF 

IIDE 

I 

J 

L 

M 

N1 

N2 

PHI 

PHID 

S 

51 

52 
Y 
Z 


A numerical constant which is the steady response of the deemphasis filter to 
a dc input voltage. 

The quotient of a component of the vector B and a corresponding component of 
the vector A. 

The peak frequency deviation. 

An increment of time T. 

A vector used to retain the sample values of cosc^(t). 

A vector used to retain the sample values of sin<£(t). 

A vector used to retain sample values of impulse response function of the pre- 
detection filter. 

A vector used to retain sample values of the impulse response function of the 
deemphasis filter. 

An integer index used in DO- Loops. 

An integer index used in DO-Loops. 

An integer index used to furnish a special ordering for vector components. 

An integer specifying the number of calculations. 

An integer specifying the number of calculations. 

An integer specifying the number of calculations to be printed out. 

A vector, the components of which are generated by the arctangent of C2. 

A vector of sample values of the derivative of PHI. 

A summation of terms involved in a convolution operation. 

A summation of terms involved in a convolution operation. 

A summation of terms involved in a convolution operation. 

A vector of sample values of the final output response function of the system. 

A vector of sample values of the partial final system output to which a single 
constant should be added to each component to give the final output vector. 
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D. 16 VARIABLES USED IN COMPUTER PROGRAM 16 

A A vector representing the convolution of cos$(t) with H^(t), the impulse response 

of the equivalent lowpass filter of a predetection filter. 1 

AB The absolute value of the derivative of PHI. 

AK The constant relating peak frequency deviation to the voltage of the modulating 

wave. 

A1 The value of an expression used repeatedly. 

B A vector representing the convolution of sin^(t) and h^(t). 

Cl A numerical constant which is the steady response of the deemphasis filter to 

a dc input voltage. 

C2 The quotient of a component of the vector B and a corresponding component of 

the vector A. 

DP The peak frequency deviation. 

DT An increment of time T. 

FLCOS A vector used to retain the sample values of cos^(t). 

FLSIN A vector used to retain the sample values of sin</>(t). 

HBF A vector used to retain sample values of the impulse response function of the 

predetection filter. 

HDE A vector used to retain sample values of the impulse response function of the 
deemphasis filter. 

I An integer index used in DO- Loops. 

J An integer index used in DO- Loops. 

L An integer index used to furnish a special ordering for vector components. 

M An integer specifying the number of calculations to be performed. 

N1 An integer specifying the number of calculations to be performed. 

N2 An integer specifying the number of calculations to be printed out. 

PHI A vector, the components of which are generated by the arctangent of C2. 

PHID A vector of sample values of the derivative of PHI. 

S A summation of terms involved in a convolution operation. 

51 A summation of terms involved in a convolution operation. 

52 A summation of terms involved in a convolution operation. 
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Y A vector of sample values of the final output response function of the system. 

Z A vector of sample values of the partial final system output to which a single 

constant should be added to each component to give the final output vector. 


f; 
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